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I attended the USACAS 2008 Conference in Chicago-Northfield which was held 28/29 June. 
Among many interesting lectures was one very special presentation. Philip Todd from 
Tigard, Oregon showed his Symbolic Geometry program Geometry Expressions. 

This is an exciting piece of software which combines dynamical geometry program with sym- 
bolic calculations in a very impressive way. 

Here are two screen shots. The first one shows the general formula of the distance of incen- 
ter and circumcenter of a triangle given by its three sides a, b and c. 


^a-jb^c 4 +c 3 -(-a-b)+c 2 -(^-a' i +3-a-b-b 2 ]+c-[a 3 -a 2 -b-a-b 2 +b 3 ] 
•/a+b+c ■Jta+b-cHa-b+cM-a+b+c) 



The second example shows the implicit form of the locus of centers of common tangents to 
two circles. 


s.4X 4 +aX 2 y 2 +4Y 4 -l2X 3 ‘a-l2XY 2 a+a 4 -a 2 s 2 +Y 2 (4a 2 -4s 2 ]4-X 2 [l3-a 2 -4s 2 ]+x[-6a 3 +4as 2 ]=0 



Find more information and download the manual (pdf) and two collections of worked exam- 
ples (pdf) at 

www.qeometrvexpressions.com 
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Dear DUG Members, 

Let me start with an excuse for being too late with DNL#70. Due to some urgent family obligations 
and my participation at USACAS 2008 in Chicago I was unable to finish this issue by end of June. 

At the other hand I am glad ot offer an issue with 48 pages content. When starting composing this 
DNL I intended to have a smaller issue but as you can see it is again full of - hopefully - interesting 
stuff. 

At USACAS 2008 I attended several interesting lectures. Among others I met Philip Todd and Steve 
Arnold. Philip demonstrated a fascinating piece of software - Geometry Expressions. ITike to invite 
you to visit his website (see Information page). Steve had a well attended excellent session on the di- 
dactical use of TI-NSpire. Have a look on his website, too and download lots of NSpire materials 
(page 3). 

Steve is one of the keynote speakers at TIME 2008. All keynote speakers and the titles of their lectures 
are given in the box below. 

You can find the abstracts of all accepted submissions for lectures and workshops for the DERIVE & 
Tl-strand and a list of all accepted lectures for the ACDCA-strand for TIME 2008. 

We have three main contributions in this DNL: Peter Schofield deals with periodic decimal numbers. 
His article reminded me how I performed calculations when I was a pupil. Benno Grabinger demon- 
strates - once more - how mathematics can be found everwhere around us. This time is a toothpaste 
tube in the focus of his interest. Finally Johann Wiesenbauer offers his Titbits 35, which are in his 
opinion one of his best Titbits. Among very sophisticated routines connected with the RSA-algorithm 
he offers a function which enables the output of calculation time for any program and/or function. (We 
had this in one of the earlier DNLs. In DNL#55 Albert Rich answered a respective request from our 
Swiss friend Rene Hugelshofer). 

I wonder if anybody can treat the simultaneous equation given on page 37 in a satisfying way with 
DERIVE. 

With my best wishes for a fine summer 
I remain 



It can happen that DNL#71 will be late, too. We will attend TIME 2008 (22 - 26 Sept.) and then have 
a 3 weeks roundtrip in South Africa. 


TIME 2008 - The Keynotes 

David Arnold : Meaningful Algebra with CAS 

Bernhard Kutzler : Technology and the Yin and Yang of Mathematics 

Education 

Nurit Zehavi : Didactical practices of computer algebra in mathematics 

education 

David Jeffrey: Debugging Computer Algebra; Debugging Mathematics. 
A Two Way Street 


Download all DNL-Derive- and Tl-files from 

http : / / www . austromath . at/ dug/ 
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Contributions: 

Please send all contributions to the Editor. 
Non-English speakers are encouraged to 
write their contributions in English to rein- 
force the international touch of the DNL. It 
must be said, though, that non-English 
articles will be warmly welcomed nonethe- 
less. Your contributions will be edited but 
not assessed. By submitting articles the 
author gives his consent for reprinting it in 
the DNL. The more contributions you will 
send, the more lively and richer in contents 
the DERIVE & CAS-77 Newsletter will be. 


Next issue: September 2008 

Deadline 15 August 2008 

Preview: Contributions waiting to be published 

Some simulations of Random Experiments, J. Bohm, AUT, Lorenz Kopp, GER 

Wonderful World of Pedal Curves, J. Bohm 

Tools for 3D-Problems, P. Luke-Rosendahl, GER 

Financial Mathematics 4, M. R. Phillips 

Hill-Encription, J. Bohm 

Farey Sequences on the Tl, M. Lesmes-Acosta, COL 
Simulating a Graphing Calculator in DERIVE, J. Bohm 
Henon & Co, J. Bohm 

Do you know this? Cabri & CAS on PC and Handheld, W. Wegscheider, AUT 

An Interesting Problem with a Triangle, Steiner Point, P. Luke-Rosendahl, GER 

Overcoming Branch & Bound by Simulation, J. Bohm, AUT 

Diophantine Polynomials, D. E. McDougall, Canada 

Graphics World, Currency Change, P. Charland, CAN 

Cubics, Quartics - interesting features, T. Koller & J. Bohm 

Logos of Companies as an Inspiration for Math Teaching 

Exciting Surfaces in the FAZ / Pierre Charland 's Graphics Gallery 

BooleanPlots.mth, P. Schofield, UK 

Old traditional examples for a CAS - what's new? J. Bohm, AUT 

Truth Tables on the Tl, M. R. Phillips 

Advanced Regression Routines for the TIs, M. R. Phillips 

Where oh Where is IT? (GPS with CAS), C. & P. Leinbach, USA 

Embroidery Patterns, H. Ludwig, GER 

Mandelbrot and Newton with DERIVE, Roman Hasek, CZ 

Snail-shells, Piotr Trebisz, GER 

A Conics-Explorer, J. Bohm, AUT 

Exercise Long Division with DERIVE 

Practise Working with times 

and others 


Impressum: 

Medieninhaber: DERIVE User Group, A-3042 Wiirmla, D'Lust 1, AUSTRIA 
Richtung: Fachzeitschrift 
Herausgeber: Mag. Josef Bohm 


The DERIVE-NEWSLETTER is the Bulle- 
tin of the DERIVE & CAS- 77 User Group. 
It is published at least four times a year 
with a contents of 40 pages minimum. The 
goals of the DNL are to enable the ex- 
change of experiences made with DERIVE, 
77-CAS and other CAS as well to create a 
group to discuss the possibilities of new 
methodical and didactical manners in 
teaching mathematics. 

Editor: Mag. Josef Bdhm 
D'Lust 1, A-3042 Wiirmla 
Austria 

Phone/FAX: 43-(0)2275/8207 
e-mail: nojo.boehm@pgv.at 
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William Pickles, Petersfield, UK will6willpwr . demon . co . uk 

Josef 

I just dug out my HP95LX running DERIVE* and can't seem to get the display running with small 
characters — is there a mail list where I might be able to get some help getting this right ? 

Regards 

William 

* I put this combination together myself about 10 years ago. 

Sorry, but I don't have any idea. I'll put your mail into our User Forum. Hopefully one of our 
members will have some advice for you, Josef. 

Wolfgang Propper, Niirnberg, Germany 

Dear Josef, 

finally I found some time to read DNL#69. 1 found a mistake: Hubert's TI-Nspire files wurf . tns 
and wurf 2 . tns (page 4) are missing in mth69.zip. 

Regards 

Wolfgang 

Another sorry. I apologize and include both mentioned files into mth70.zip. Thanks for the 
note, Josef. 


Reanimating a TI-89 

Renate Wronski, a colleague from Styria wrote that her TI-89 didn't show any sign of life af- 
ter having a break of some months because of her sabbatical. As she didn't receive any an- 
swer from Tl I offered to send her one of my Tl-89s. But then she wrote back: 

Renate Wronski, Graz, Styria 

I tried to reach the Tl-service by phone and had success. The following procedure brought my TI-89 to 
life again: 


Remove one of the four batteries. Keep the APPS-key pressed while setting in this battery 
again. The device will 'wake up'. Then remove one battery again and and set it in without 
pressing any key. You should see now the black bar on the screen which is the sign for 
loading the OS. The you can recognize the well known TI-89 desktop. 

This is a way to reset your calculator. I hope that there will be no more troubles in the future. 

Best regards and many thanks for your offer, 

Renate Wronski 


For our TI-Nspire Users: 

At the occasion of USACAS 2008 Conference in Chicago-Northfield I had the occasion to 
meet Stephen Arnold from Australia who is one of the TIME 2008 keynote speakers. Steve 
has a great website with a bundle of TI-NspireCAS papers: 


http : //compasstech . com.au . 
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Piterr, Poland playmakerpit 0 poczta . onet . pi 

Hello 

how calculate this differential equation: 


4/-Vj-l = 0. 

Thanks, Piterr 

Dear Piterr, 

This is the result of the DERIVE attempt - according to the Online Help 
(see under Second order Ordinary Differential Equations how to interpret the general solu- 
tion): 


DERIVE Online Help: 

AUTONOMOUS_CONSERVATIVE(q, x, y, xO, yO, vO) simplifies to an im- 
plicit algebraic solution of an equation of the form y"=q(y) having initial condi- 
tions y=yO and y'=vO at x=xO. Equations of this form are autonomous be- 
cause the variable x does not occur and conservative because y' does not 
occur. If instead of the initial condition y'=vO you have a second boundary 
condition y=y2 at x=x2, substitute x2 for x and y2 for y in the solution and 
solve for vO, which you can then eliminate in favor of x2 and y2. 


#1: AUTONOMOU S_CON S E RVATIVE 

l 4-7y 

2 

4* (2 *v0 - B-^yO). |v0| -IF(vO = 0, 1, SIGN(vO)) 

#2: x = + 

3 


2 2 

4*7(7y + vO - 7y0)-(7y - 2-(v0 - 7yO))-IF(vO = 0, 1, SIGN(vO)) 


The Tl-Voyage 200 has a built-in DE-Solver, which returns the general solution together with 
two constants @1 and @2. 


|Fi 7 »ir|Y FIT Y F3t Y FHt Y FE Y FbT Y 'I 
|t |fllgebra|Calc|0ther|PrgnI0|ciean Up| f 

■ deSol 1 

■ solce| 
y C-Jy ■ 

ue(4 ■ y 1 1 ■ Jy - 1 = 0 * X , y ) 

4-J35Tm-(J5-2-e2) 

3 =x + ei 

f 4-Jfy + (22 (fy -2(22) _ ] 

L 3 -x + Gl.yJ 

, 9 x 2 + IS ■ ei x + 9 ■ ei 2 - 64fc 

-3e2j- 16 

i 

DC 

RAD EXACT FUNC 1H/30 


Answer from Piterr: 

Tha nks 

Function AUTONOMOUS_CONSERVATIVE(q, x, y, xO, yO, vO) must have begin condition. My 
example don't have begin condition. 

Hm. . . Could I use this function ? I review Derive help and I can't match function for this example. 
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Abstracts of the ACCEPTED proposals for TIME 2008 - 
DERIVE & TI-CAS - Conference 

1 Revisiting surprising results with CAS calculators 

(Lecture 25 min) 

Gilles Picard and Michel Beaudin, Ecole de technologic superieure, Montreal, Canada 

We teach a variety of math topics (review of College Algebra, Calculus, Differential Equations, ap- 
plied probability and statistics) in a Technical Engineering School. The Voyage 200 (or TI 89 Tita- 
nium) is mandatory for all new full-time students. We make use of this calculator on a regular basis, 
for exploring with students, in the classroom, all the classical curriculum in mathematics. 

In July 2006, in Dresden, using the Voyage 200 handheld calculator, equipped with the operating sys- 
tem OS 3.10, we showed some surprising results given by the device. We showed examples where the 
CAS system gave unsatisfactory or strange results or where it couldn't perform some commands. 
These problems were often related to the way the CAS system would simplify (or not) some expres- 
sions in some intermediate step of the calculations. We made some suggestions in order to fix the en- 
countered bugs. We are happy to see that in Nspire CAS, OS 1.3, some of these bugs have been fixed 
and, in this talk, we will use Nspire CAS (both handheld and software) to show it and to see which 
bugs still remain. After revisiting these old examples, we will show some new problems encountered 
with our students while using the Voyage 200 CAS engine. These surprising results are still not re- 
solved as of Version 1.3.2437 (2008-01-08) of the Nspire CAS calculator. 

2 Functions, Programs and Libraries with TI-Nspire 

(Workshop 90 min) 

Josef Boehm, ACDCA, A 3042 Wuermla, Austria 

TI-NspireCAS version 1.3 offers a very comfortable program editor. In this workshop we will focus 
on the difference between functions and programs by introductory examples. The second issue will be 
demonstrating how to create and to work with libraries which can be used very similar to utility pack- 
ages in other computer algebra systems. TI-Nspire units will be provided. 

3 Exploring Zeros of Complex Functions Graphically 

(Lecture 25 min) 

Giora Mann and Nurit Zehavi, Weizmann Institute of Science, Israel 

Implicit plotting - one of the tools available in CAS - makes the exploration of zeros of any complex 
function a very simple procedure, as we shall demonstrate. Historically, the constraint on solving 
complex equations was that the graph of a complex function is 4D. Our idea was that both loci can be 
plotted implicitly, and the zeros are the intersection points of the two curves. However, not all the 
points that look like intersection points are actually intersection points. As a measure of control we 
introduce contour maps of |f(z)| 2 to determine whether a ’suspected’ intersection point is indeed an 
intersection point, which is a zero of the function. 
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Students who are familiar with analysis of 2-variable functions can use their knowledge for finding the 
minimum points of |f(z)| 2 as the intersection points of Re(|f(z)| 2 ) and Im(|f(z)| 2 ). There are two draw- 
backs in this approach: (a) the students need to understand partial derivatives and (b) zero partial de- 
rivatives is only a necessary condition; it means that we get sometimes as many saddle points as 
minimum points. 

The idea of zero imaginary and real parts of a function at a certain point, making this point a zero of 
the function, is a basic idea of complex analysis, which makes it independent of 2-variable analysis. 
This procedure is general and can be integrated into an introductory course of complex analysis, re- 
quiring only basic knowledge of 2-variable real functions and complex numbers. 

4 Real-life applications of ODEs for undergraduates 

(Lecture 25 min) 

YuHe Yuan, Steve Joubert, Ying Gai., Dept, of Mathematics and Statistics, Tshwane Univer- 
sity of Technology (TUT), Pretoria, South Africa 

This study introduces real-life mathematical theories and models of international relationships suitable 
for undergraduate ordinary differential equations, by investigating conflicts between different nations 
or alliances. Based on the work of Richardson, systems of differential equations are constructed. The 
solutions and the stability of systems of ODEs are observed, with the aid of mathematical softwares 
such as Derive, Mathematica and Scientific Workplace. One of the most interesting tasks is to analyse 
the coefficients in the constructed models. In our opinion, the model first constructed by Richardson is 
an excellent application of ODEs (ordinary differential equations) and is useful for practice for learn- 
ing ODEs. One would expect this kind of model to be added to the material in textbooks as a typical 
example. 

5 Can CAS be trusted? 

(Lecture 25 min) 

Stephan V. Joubert and Temple H. Fay, Department of Mathematics and Statistics; TUT, Pre- 
toria, South Africa 

Most computer algebra systems (CAS) have built-in ordinary differential equation (ODE) solvers, but 
the accuracy of the solutions produced is not always obvious. Various ways of estimating the accuracy 
of ODE solvers are discussed here, extending work presented at the ’’Remarkable Delta 2003” confer- 
ence in Queenstown, New Zealand and the ’’TIME 2004 Conference” in Montreal, Canada. Our meth- 
ods are easy enough for undergraduates to implement because the needed mathematics is accessible to 
them. Many students (and their teachers) have an in-depth knowledge of how to check the accuracy of 
numerical routines, but many trust them blindly. On the other hand, testing the accuracy of a routine 
takes more time than just running the routine to produce a solution and this is another reason for tak- 
ing a solution at face value. Such blind trust could have negative connotations if carried through to 
industry and elsewhere after the student graduates. We cite an example of how experienced mathe- 
matical scientists (academics) have fallen into the trap of assuming numerical solutions to be correct. 
There already exist a number of routines to test the accuracy of ODE solvers, some of them time in- 
tensive, and some not. The routines introduced here add to this collection of routines and one of them 
substantially reduces the calculation time of an existing routine. We extend our results for initial value 
problems to boundary value problems in ODE. 
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6 Data Acquisition and mathematical modelling - A case study 

(Workshop 90 min) 

Anna C.M. Bekker, Stephan V. Joubert and Temple H. Fay, TUT, Pretoria, South Africa. 

A mathematical model is derived of a motor-vehicle tyre tread surface striking a speed bump on a 
stretch of otherwise smooth horizontal road. A mathematical idea is outlined to derive the model and a 
simple experiment, involving a typical vinyl record player needle, is described to measure the coeffi- 
cient of damping of the tangential vibration of the rubber tyre. 

7 A Discreet Compartmental Model for Lead Metabolism in the Human Body 

(Lecture 25 min) 

Charlotta E Coetzee, Stephan V Joubert and Frederika E Steyn, TUT, Pretoria, South Africa 

A real-life example of a mathematical technique, employing a so-called transfer matrix, is developed 
for the compartmental analysis of lead metabolism in the human body. The technique is demonstrated 
with the aid of Bert's four-compartment biokinetic model. The results produced by this time-discrete 
approach correspond almost perfectly with those of a continuous-time method. The powerful calcula- 
tion tools of the Computer Algebra System (CAS) "Scientific Workplace" are employed to illustrate 
the results using tables and graphs. 

8 Expanding Student Perspectives: A Workshop on Forensic Applications of 
Mathematics 

(Workshop 90 min) 

Patricia Leinbach and Carl Leinbach, Adams County Pennsylvania Coroners Office and Get- 
tysburg College, PA, USA 

In this workshop you will be a member of a crime scene investigative team conducting a forensic in- 
vestigation of the scene. Your job will not be to solve the crime, but to determine the manner and 
cause of death of an individual found at the scene. You will gather and analyze "evidence" gathered at 
the scene. (In fact, because this is a workshop on using such investigations in the classroom, you may 
be asked to create some "evidence" related to the description you are given.) The workshop will con- 
clude with your team presenting their results and analysis to their colleagues at the workshop. 

9 Using Learning Objects with TI-Nspire CAS 

(Workshop 90 min) 

Wade Ellis, Jr., West Valley College, San Jose, CA, 95130, USA 

In this workshop, participants will work with a variety of learning objects developed for TI-Nspire 
CAS. They will also examine and discuss activities (lessons) that incorporate these learning objects. 
A learning object is a TI-Nspire document that allows students to act on a mathematical object, ob- 
serve the consequences of their actions, and then reflect on the mathematical meaning of those conse- 
quences. These learning objects are intended to be use as investigative tools by students to generate 
and enhance mathematical understanding. These objects provide a platform for the development of 
investigative activities as well as problem-solving activities. 
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10 A Trinomial Factoring Investigation with Pre-service teachers 

(Lecture 50 min) 

Michael Meagher, Michael Todd Edwards, Asli Ozgun-Koca, Brooklyn College - CUNY, Mi- 
ami University, Wayne State University, USA 

Pre service teachers tend to accept the secondary mathematics curriculum as a static set and tend not to 
question the inclusion or exclusion of certain topics. Factoring of trinomials is a standard topic in 
school mathematics curricula worldwide and much class time is devoted to this topic. Building on 
work from Usiskin we present an activity, using the CAS and spreadsheet capabilities of the TI- 
Nspire, which engages students in establishing what percentage of trinomials with integer co-efficients 
are factorable. The task, which involves algebra, probability and statistics was then used as a vehicle 
to discuss with students the inclusion of specific topics in the school curriculum and how such deci- 
sions are made. The presentation will engage participants in the activity itself and will report the re- 
sults of the students' discussion of the curriculum. 

11 Separatrices 

(Lecture 25 min) 

Phindele M. Skhosana, Temple H. Fay and Stephan V. Joubert, TUT, Pretoria, South Africa 

In this presentation we examine two by two first order systems of ordinary differential equations and 
show how to determine phase plane portraits and identify separatrices when there is a saddle point. In 
order to do so we describe how to use a computer algebra system to generate trajectories from contour 
plots when possible and from numerical investigations. In many cases we can determine the equation 
of the separatrix. Generating a phase plane portrait is useful, for at a glance one can observe what ini- 
tial values give rise to bounded solutions, periodic solutions and other important features. It also per- 
mits the instructor to concentrate on the qualitative aspects of the model under investigation rather 
than the calculational difficulties associated with finding solutions. 

12 Modelling Cha Cha dance in using the "function"-tools within the Cabri 2 Plus 
or the TI NSpire environment 

(Lecture 50 min) 

Jean-Jacques Dahan, IREM of Toulouse, France 

We will show how to model Cha Cha dance in programming the movement of two points (modelling 
the two feet of the dancer). We will use the special features of Cabri 2 plus and TI Nspire to lead these 
points with functions defined on different intervals. We will compare the two different approaches: 

The one with Cabri where it is possible to superimpose the graphic and the algebraic frames on the 
same page. 

The one with TI Nspire where it is possible to display these different frames on different pages in the 
same screen 

13 TI-Nspire Learning Technology 

(Lecture 50 min) 

Gosia Brothers, Texas Instruments, Dallas, TX, USA 
Come to our presentation to learn more about TI-Nspire Technology. 
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14 Improving Algebraic Environment 

(Workshop 90 min) 

Gosia Brothers, Texas Instruments, Dallas, TX, USA 

Come to our session to give us your feedback on TI-Nspire Algebraic Environment and discuss your 
wishes for the future versions of products from Texas Instruments. 

15 Improving Graphs and Geometry 

(Workshop 90 min) 

John Good, TI-Team 

Come to our session to give us your feedback on TI-Nspire Graphs and Geometry and discuss your 
wishes for the future versions of products from Texas Instruments. 

16 Some Important Functionalities of a CAS when Teaching Mathematics 
to Future Engineers 

(Lecture 25 min) 

Michel Beaudin, ETS, Montreal, Canada 

We are teaching mathematics at Ecole de technologie superieure (ETS), an engineering school in 
Montreal, Canada where every student has a Voyage 200 calculator on his desk and has access to 
computer labs where CAS like Derive, Maple and also Matlab program are installed. In Vienna 
(Visit-me 2002), we showed many examples of how Derive 5 and the TI-92 Plus were used when 
teaching to engineering students. For the past 4 years, we used the Voyage 200 and Derive 6.10 for 
teaching single and multiple variable calculus, linear algebra, differential equations, complex analysis. 
The talk will show examples of the importance of 2D implicit plotting, 3D plotting and Odes plotting 
when teaching to future engineers. These features are not yet implemented into Nspire and we hope 
that it will be on board soon. If we agree that we have to make a move from Derive to Nspire, we 
don't agree to leave higher mathematics subjects to the competitors. 

17 Nanopowder Production in a Plasmachemical Reactor: 

Computational Fluid Dynamics Modelling and Simulation. 

(Lecture 25 min) 

Phethedi J. Kekana; Andrei V. Kolesnikov; Stephen V. Joubert., TUT, Pretoria, South Africa 

In this present study the development of a more realistic two-dimensional mathematical model capable 
of predicting the main aerosol phenomena such as Ti02 nanoparticle formation, growth and deposi- 
tion in a high temperature plasmachemical reactor and the simulation results are presented. The Ti02 
nanoparticles were produced in a tubular reactor from the oxidation of TiC14 vapor in an Argon at- 
mosphere. The developed mathematical model consists of mass, momentum and energy conservation 
equations. The particle dynamics processes include particle formation by nucleation, growth by con- 
densation and coagulation as well as the loss of product particle by deposition on the wall of the plas- 
machemical reactor. The aim of this model was to predict both the axial and radial profiles of the flow 
velocity, temperature and concentrations profiles of TiC14, 02, Ti02 and C12 and most importantly 
the evolution of the particle size distribution of Ti02 nanoparticle and Ti02 nanoparticle deposition to 
the reactor wall by numerical integration of stiff nonlinear partial differential equations 
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(PDE) in the FLUENT CFD program incorporating the Method of Moment technique. In this model 
the particle size distribution was approximated by the Gaussian lognormal function. The simulation 
results obtained from a two-dimensional model provides us with useful information on the influence of 
operational conditions (i.e. gas components flow rates, initial temperature and concentrations) and the 
reactor configurations on the evolution of Ti02 particle size distribution and the deposition flux at 
reasonable computational time and memory. The performance of the detailed two-dimensional model 
was validated by comparing its predicted results with the experimental and/or numerical data already 
published in the literature. Good correlations between the predicted and experimental results were 
achieved. 

18 Using Science as a tool for learning mathematics 

(Lecture 25 min) 

Hildegard Urban, Vienna, Austria 

Mathematics, as the language of numbers, is an important tool in science classes, but science is not 
generally considered as a tool for teaching mathematics. This article presents examples incorporating 
science concepts and problem solving in math classes using a motion detector (Calculator Based 
Ranger, CBR) and technology from Texas Instruments (TI-Nspire-CAS-handheld or. TI-Nspire-CAS- 
computer software). Real world data collection tools and Nspire introduce students to many fascinat- 
ing concepts in mathematics and give them interactive ways to visualize relationships and patterns and 
enhance critical thinking. The author is investigating the mathematical and pedagogical potential of 
using TI technology (Graphical calculator, Voyage, Nspire) in combination with Vernier sensors and 
probes as devices to collect various kinds of data and of using the software to serve as a powerful 
analysis tool, helping students to build mathematical models. Experiences have been made in grade 9 
to 11 (15- to 17-year old students) are reported. The use of technology seems to effectively enhance 
students' learning. Students are actively engaged in learning as they make predictions, take measure- 
ments, analyze their data and make decisions about presenting their work. They are challenged to dis- 
play their individual talents and mathematical abilities in real world problem solving situations. 

19 The Many Dimensions of Decision Making A Process for Making Decisions 

in a Complex Environment 

(Lecture 50 min) 

Carl Leinbach, Professor Emeritus, Gettysburg College, PA, USA 

In today's society no sound decision is made simply on the basis of one issue alone. Furthermore, 
reasonable people, looking at the situation may make entirely different decisions. There is no "one 
size fits all." The reason is that different people and different groups have different priorities and de- 
sires. The question we will consider is, how a decision maker can bring together these differing priori- 
ties into a group consensus. This presentation will present a process based on Thomas Saaty's Ana- 
lytical Hierarchical Process for making complex decisions. It will begin with a brief discussion of the 
underlying mathematical assumptions of the process and apply the technique to the problem of choos- 
ing an afternoon tour at the TIME 2008 Conference. If time permits, the audience will participate in 
constructing a group priority for an issue of current interest. 
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20 Teaching Mathematics to Engineering Students: To Use or 
Not To Use TI-Nspire CAS 

(Lecture 25 min) 

Michel Beaudin and Gilles Picard, ETS, Montreal, Canada 

We have been using the Voyage 200 to teach a variety of math topics in a Technical Engineering 
School (ETS in Montreal, Canada). We have to admit that the Voyage 200 is doing a very good job 
but, for some problems, a much faster processor would be helpful and would allow more powerful 
graphic options. This fast processor is exactly what we now have in the TI-Nspire CAS. So, when you 
need to solve some "heavy problems”, say complicated equations, polynomial systems of equations, or 
when you want to define special functions using a definite integral, the Voyage 200 processor cannot 
compete with the TI-Nspire CAS calculator. The talk will give some examples of this, using Nspire 
CAS and comparing with Voyage 200 results. 

Since we are teaching mathematics to future engineers, we also often need implicit 2D plots, 3D plots 
and differential equations plotting (along with RK and Euler numerical methods). We have all of this 
in the Voyage 200 but it can take a lot of time to get some results or see graphs on the screen and, in 
3D plotting, the results are limited mainly due to a lack of processor speed. Adding these features to 
Nspire CAS, and colour for the PC version, would be important and, as far as we are concerned, this is 
a must. If this is not done, 3D plots will continue to be done, by our colleagues and us, using Derive 
or Maple software. In order to make the move to Nspire CAS software, we absolutely need a much 
more "university level package". 

21 Modelling of the telegraph equations in transmission lines 

(Lecture 25 min) 

Corrie Lock, Proff JC Greeff, SV Joubert UJ, TUT, TUT, South Africa 

It would be difficult to imagine a world without communication systems. A plethora of guided fixed 
line telephones as well as a multitude of unguided systems to serve cellular phones are evident in our 
surrounding world. In order to optimise guided communication systems, it is necessary to determine or 
project power and signal losses in the system, since all systems have such losses. To determine these 
losses and eventually ensure a maximum output, it is necessary to formulate some kind of equation 
with which to calculate these losses. A mathematical derivation for the telegraph equation in terms of 
voltage and current for a section of a transmission line will be investigated. 

In literature for engineers consulted, the formulae for voltage and current involved in the telegraphic 
equations are not explicitly and analytically derived, leaving a theoretical gap seldom crossed by stu- 
dents in Electrical Engineering. 

The main aim is to address this theoretical gap, and derive from basic principles, the equations for 
telegraphic transmission in a guided system and secondarily to illustrate the applications thereof to 
real-world problems using a suitable computer algebra system in this case, Derive. 
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22 Comparison of an Analytical Method and Matlab to Model Electromagnetic 
Distribution in a Trough 

(Lecture 25 min) 

JJ Bruyns, J.C. Greeff, S.V. Joubert, UJ, TUT, TUT, South Africa 

In designing devices that discharge electrostatic voltage it is important to know the radiation pattern. 
To interpret and visualize the result is just as important as being able to obtain a correct analytical 
result. This lecture discusses the traditional analytical method to determine the radiation pattern in a 
closed trough, as well as the use of a computer algebra system (CAS) to determine the radiation pat- 
tern. MATLAB is used as a computer algebra system to solve and visually display the electrostatic 
distribution in a trough. 

The results of the analytical and numerical solutions are compared to determine the accuracy of the 
CAS solution. The use of the CAS system and its educational advantage is explained. 

23 Teaching differential equations and its application - Using Derive 6 as a PeCAS 

(Lecture 25 min) 

Jose Luis Galan, Gabriel Aguilera et al, University of Malaga, Malaga, Spain 

In this talk we will describe the file DIFFERENTIAL EQUATIONS.MTH, created in DERIVE 6 in 
order to be used in mathematical subjects which deal with differential equations, aimed at Engineering 
students. Such file contains a series of programs which permit to solve differential equations problems. 

The programs contained in the file can be grouped within the following blocks: 

First order differential equations: separable equations and equations reducible to them, homogeneous 
equations and equations reducible to them, exact differential equations and equations reducible to 
them (integrating factor technique), linear equations, the Bernoulli equation, the Riccati equation. 

• First order differential equations and nth degree in y\ 

• Generic programs to solve first order differential equations. 

• Cauchy problems for first order differential equations. 

• Higher orders differential equations. 

• Cauchy problems for higher orders differential equations. 

• Applications of differential equations. 

We will also show in the talk some examples of applications that have been carried out with our stu- 
dents of Telecommunication Engineering. 

The programs have been developed using the Display function in order to be used as didactical tools 
with explications of what the programs do step by step, using DERIVE 6 as a Pedagogical CAS (Pe- 
CAS) or as a white-box CAS. 

Finally, we include the conclusions obtained after using this file with our students and also some fu- 
ture work on this subject. 
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24 DERIVE 6 as a pedagogical CAS: using the slide bar utility and the 
DISPLAY function in programming 

(Workshop 90 min) 

Jose Luis Galan, Gabriel Aguilera et al, University of Malaga, Malaga, Spain 

In this workshop we will use some examples of programming with DERIVE 6 using the display func- 
tion in order to use this software as a pedagogical CAS (PeCAS). We have developed this kind of 
workshop with our students of Technical Telecommunication Engineering. 

The main innovative aspect of this way of teaching is that students have an active role. Specifically, 
they have to elaborate, by themselves, utility files to solve typical problems for different subjects. In 
our case, this fact implies that students need to deal with programming in DERIVE 6, understand the 
subject and know how to solve typical problems. 

This workshop will consist of two different parts: 

1 . Drawing classical curves using the slide bar utility. 

In this first part, we will draw curves such as segment, circumference, ellipse, lemniscate, astroid, 
cardioid, catenary, cycloid, cissoid, folium, eight curve, limacon, rhodonea curves, hypocycloid, tri- 
sectrix, tractrix, spiral, . . . 

Changing, with the slide bar utility, the parameter(s) of a curve, we can study properly its properties. 

2. Programming and computing line integrals using the display function 

The second part of the workshop will consist of the development of different programs in Derive 6 in 
order to solve any line integral considering both, computing the integral by means of the definition 
and/or using an appropriate theorem to compute it properly. 

We will develop the programs using the display function which will allow us to use DERIVE 6 as a 
PeCAS. 

We will provide some material to make the workshop easier to follow. 

25 New Kids on the Block - Some First Experiences with Recent Alternatives to DERIVE 

(Lecture 25 min) 

Karsten Schmidt, Schmalkalden University of Applied Sciences, Germany 

Although Texas Instruments finally discontinued DERIVE last summer, numerous DERIVE users are 
still working with this popular Computer Algebra System for a variety of reasons. But in addition to 
the well-known DERIVE competitors - Mathematica, Maple, and MuPAD - there are now new alter- 
natives in the Computer Algebra System market, such as Nspire and WIRIS. 

In this presentation we reflect on possible difficulties in the transition process, from a more organiza- 
tional point of view (Is the licensing system appropriate? What is the price tag? How easy is the instal- 
lation?), as well as from a didactic point of view (How intuitive is the software for inexperienced us- 
ers? How much class time is required before the respective Computer Algebra System can be used 
effectively?). 

These considerations are made from the background of a first-year university course in linear algebra 
which was transformed step by step over a period of ten years from a traditional "blackboard and 
transparencies” teaching approach (in a lecture hall standing in front of up to 70 students) into an in- 
teractive teaching approach using DERIVE (in a PC lab with 20 PCs and no more than 40 students). 
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We will also be looking at the question as to whether the new kids on the block have something inter- 
esting to offer in this context which is not available in DERIVE. 

26 Introducing a computer algebra system in mathematics education - empirical evidence 
from Thuringia (Germany) 

(Lecture 50 min) 

Dr. Wolfgang Moldenhauer, Lehrplanentwicklung und Medien, Heinrich-Heine-Allee 2-4, 
99438 Bad Berka, Germany 

This lecture reports on the effects the use of a pocket calculator-based computer algebra system (CAS) 
has on the performance in mathematics of grade 11 students in Thuringia. A project started at 8 of 
about one hundred upper secondary schools in the federal state of Thuringia in 1999; 3 years later the 
former restrictions on the use of technology in math education were lifted. In 2004, more than a quar- 
ter of all Thuringian upper secondary schools used CAS in math classes. Beginning in 2000, a test was 
carried out each year to compare the performance of CAS and non-CAS students (from different con- 
trol schools). More than 12000 students were tested. In 70% of the cases CAS students performed 
better than, and in the remaining 30% they performed as well as, non-CAS students. There is evidence 
that students in advanced courses benefit more from using CAS than students in basic courses. 

27 A novel method of interpolation and extrapolation of functions by a 
linear initial value problem 

(Lecture 50 min) 

Michael Shatalov, Igor Fedotov and Stephan V. Joubert, CSIR and TUT, South Africa 

The classical approach to function approximation is based on a particular choice of functions, for ex- 
ample polynomial, rational, exponential functions or Fourier series. The main advantage of these 
methods is to obtain the approximation expressions in a closed form. There are several disadvantages 
to the classical approach. For example, polynomial interpolation may seldom be used for the purposes 
of extrapolation due to the fast divergence of higher order polynomials outside of the interpolation 
interval. The main disadvantage of a Fourier series approximation is that it is not applicable to non- 
periodic functions and hence, could not be used for extrapolation purposes, et cetera. The method we 
propose allows us to approximate functions by means of linear combination of polynomials, trigono- 
metric and exponential functions, products of polynomials and exponents, polynomials and periodic 
functions, periodic functions and exponents, and polynomials, exponents and periodic functions et 
cetera. It is well suited for the purposes of interpolation and extrapolation of physical and chemical 
processes, which are described in terms of systems of linearized ordinary differential equations 
(ODE). The main idea of the proposed method is the approximation of a function on a fixed interval 
by means of linear ODE with unknown constant coefficients. Initial values of the problem are also 
considered as unknowns. The goal function is formulated as a positive definite function with non- 
negative weight function. Unknown coefficients and initial conditions are defined by means of mini- 
mization of the goal function. Examples of practical approximation of functions are considered and 
compared with available commercial algorithms of interpolation, extrapolation and smoothing. The 
methods we discuss will be readily understood by undergraduate students that have been taught ODE 
using DERIVE or some other CAS. 
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28 A CAS Approach to Understanding from Beginning Algebra to 
Advanced Calculus and Abstract Algebra 

(Lecture 50 min) 

William C. Bauldry and Wade Ellis, Appalachian State University (NC) and 
West Valley College (CA), USA 

The presenters will use TI-Nspire CAS to demonstrate ways to foster and enhance student understand- 
ing of mathematics in courses from secondary school Beginning Algebra to university Linear Algebra, 
Advanced Calculus, and Abstract Algebra courses. The presentation will involve the Ac- 
tion/Consequence/Re flection Principle developed by Tom Dick and Gail Burrill whereby students act 
on mathematical objects, transparently observe the consequences of their actions, and then reflect on 
the mathematical meaning of those consequences. 

29 Heat transfer in a one dimensional domain of variable cross-sections 

(Lecture 25 min) 

RS Lebelo, I Fedotov, M Shatalov and HM Djouosseu Tenkam, TUT, Pretoria, South Africa 

In this paper the method of approximating solutions of partial differential equations with variable co- 
efficients is studied. This is done by considering heat flow through a one dimensional domain model, 
with variable cross-sections of N sections. The heat transfer process is described by heat equation. 
This study is based on finding eigenvalues using software such as Derive and Mathcad. The corre- 
sponding eigenfunctions automatically satisfy the boundary conditions at the endpoints and boundary 
conditions of the first kind at the endpoints will be considered for the conic section. The authors show 
how a student can solve different eigenvalues and eigenfunctions using above mentioned softwares. 
This is a central point of finding the analytical solution of partial differential equations. 

30 Dipstick Readings 

(Lecture 25 min) 

Rambane D.T., TUT, Pretoria, South Africa 

In garages fuel (petrol or diesel) is often stored in cyclindrical tanks underground. To measure the 
volume of fuel in the tank a deepstick is used. We investigate the mathematics involved and how De- 
rive can be used. 

31 A probabilistic approach to function approximation 

(Lecture 25 min) 

PH Kloppers, TH Fay, SV Joubert, TUT, South Africa 

In this talk we will investigate a new and novel approach to function approximation for functions de- 
fined over a bounded real interval which, without loss of generality is assume to be the unit interval 
[0,1]. This approach is based upon an idea by J. Kolibal and C. Saltiel. We will give a derivation of 
their technique and show how to interpret it statistically using the Gaussian distribution. Since Bern- 
stein polynomials lie at the foundation of the Kolibal- Saltiel technique which, they call the Bernstein 
function approximation, we will call our technique the Gauss-Bernstein approximation technique. 

These ideas are quite general and have wide applicability in function approximation, high frequency 
filtering, data interpolation and multidimensional data regularization. 
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This technique can be taught at the undergraduate level because the use of a computer algebra system 
such as Mathematica is essential for the purposes of visualization and easing any algebraic, technical 
barriers that students may encounter. 

32 Deep Learning and Fun in First Year using Maple 

(Lecture 50 min) 

Bill Blyth and Alexandra Labovic, RMIT University, Melbourne, Australia 

Mathematics educators and their students need to embrace technology. Studies have shown that some 
students have some negative attitudes towards using computer packages and programming. Our objec- 
tive is for our graduates to be expert users of the powerful CAS (computer algebra system), Maple. 
Thus we want to ensure that students' initial experience is positive - particularly in the first semester 
when any negative attitudes need to be overcome. 

In the first semester of a traditional calculus course, the weekly Maple lab sessions are not used to 
directly support the lectures ... nearly half of the work is at school level! The purpose is for the stu- 
dents to enjoy the experience of using Maple. 

Student work in groups of size 2 to 4. After Maple introductions, they complete an Introduction to 
Animation session and then choose an extended animation project from a list of five problems. They 
have to demonstrate their animations in the lab for assessment. Students enjoy the animation project. 

Following the animation projects, Spot the Curve uses plots and animations to understand horizontal 
and vertical translation of curves: students identify the translation used and appreciate automatic mark- 
ing within Maple. Student feedback has been very positive. 

Our trapezoidal rule assignment is now more fun: it's disguised as a Fish Pond (a trout farm). Trape- 
zoidal rule is used to approximate the cross-sectional area (and hence the number of trout) - the stu- 
dents download a template for individualized fish ponds, with automatic marking. These projects are 
enjoyable deep learning activities. 

33 Remarks on Duffing 1 s Equation 

(Lecture 50 min) 

TH Fay, TUT, South Africa 

We discuss some of the interesting features of the forced Duffing equation 

d/dt(dx/dt) + k*dx/dt-a*x + x A 3=F*cos(b*t) 

through numerical investigations using a computer algebra system (in our case Mathematica version 
5.0). We discuss the unforced damped and undamped equation briefly, and concentrate on solving 
numerically the differential equation for a variety of values of the parameter F holding all other pa- 
rameters fixed and generally holding the initial conditions for x(0) and x(0) to be the "at rest" condi- 
tions of zero. In doing so, many striking and fascinating trajectories representing interesting motions 
and other phenomena can be discovered including: stability, periodic solutions (both harmonic and 
subharmonic), almost periodic solutions, and aperiodic solutions. In particular, chaos is often claimed 
to be evident in the trajectories and solutions of this Duffing equation and it is the purpose of this arti- 
cle to elaborate on this. These studies naturally give rise to computer laboratory problems suitable for 
student research and small group projects. Numerical investigations should go hand-in-hand with theo- 
retical studies as the one cross fertilizes the other. As an Addendum, a list of student research prob- 
lems is attached. 
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Accepted Papers for the ACDCA strand of TIME 2008 

Introduction to TI-Nspire CAS, Bernhard Kutzler , Workshop 90 min 
Exercising Control: Didactical Influences, Kathleen Pineau , Lecture 50 min 
CAS and calculation competence of students, Rainer Heinrich , Lecture 50 min 
Basic Skills and CAS, Josef Bohm, Lecture 50 min 

Linking geometry, algebra and calculus with GeoGebra, Josef Bohm, Lecture 25 min 

Modelling Cha Cha dance with Cabri 3D, Jean-Jacques Dahan, Workshop 90 min 
(was moved to the DERIVE & Tl-conference in order to accomplish Jean-Jacques lecture on Model- 
ling Cha Cha Cha with CabriII+ and NSpire) 

Spreadsheets and Interactive White Boards in the Primary Mathematics Classroom, 

Philip Oostenbroek, Br Adrian Story, Dr Anne Williams (Lecture 50 min) 

Why Proof in Dynamic Geometry?, Michael de Villiers, Lecture 25 min 

Argumentation Schemes and the Use of Sketchpad, Angel Homero, Lecture 25 min 

Neither a Tractor, nor a Matrix but a Tractrix!, Susan Steyn, Lecture 25 min 

Sustainability of mathematics education by using technology demonstrated with the 
topic of exponential growth, Helmut Heugl, Lecture 50 min 

Tasks in Calculus: Results of a 9-Year Evolution, 

Genevive Savard and Kathleen Pineau, Lecture 25 min 

Game proramming / a "future” method to teach mathematics, Sofia Backstrom & Marie Rudenstam, 
Lecture 50 min 

CAS-exercises during the Central Examination in North Rhine-Westphalia (Germany), 

Dirk Warthmann, Lecture 25 min 

Roots of transcendental algebraic equations: A method of bracketing roots and selecting initial estima- 
tions, J.N. Mwambakana, M. Shatalov, and I. Fedotov, Lecture 25 min 

Application of eigenfunction orthogonalities to vibration problems, H.M Djouosseu Tenkam, I. Fe- 
dotov, M. Shatalov, Lecture 25 min 

A School-Oriented Review of Computer Algebra Systems for Solving Equations and Simplifications. 
Issues of Domain, Eno Tonisson, Lecture 50 min 

Defining a stability boundary for three species competition models, Quay van der Hoff, Johanna C. 
Greeff, Lecture 25 min 

Technology - can it be trusted?, Johanna C Greeff, Stephan V Joubert, Lecture 25 min 

Identification of Dynamical Systems Parameters from Experimental Data using Numerical Methods, 
Pete N.A and Fedotov I, Lecture 25 min 

Numerical Computation of Special Functions with Application to Physics, Motsepe KA; Fedotov I; 
Shatlov M., Lecture 25 min 

An error analysis of the numerical method of lines, Judith N.M. Bidie, Temple H. Fay, 

Stephan V. Joubert, Lecture 25 min 

Using Matlab for teaching Mathematical Modelling, Ansie Harding, Lecture 25 min 
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For example, Simplify DtoQ(“1 .2’3456”) to obtain: 

5555 

I have also made use of a simple flag variable (default value = 0) so that, when 
flag = 1 , the procedure attempts to display the arithmetic structure of the calculation. 

2 3456 

For example, Simplify DtoQ(“1.2’3456”,1) to obtain: 1 + + 

10 99990 


2. Quotients to Decimals in Derive 6 

To form a decimal expansion of a fraction the following procedure uses the standard method of divid- 
ing of the denominator into a sequence of numerators (scaled up by base 10 each time). If (at some 
stage) the remainder becomes zero the decimal terminates, if a non-zero remainder is repeated then the 
decimal becomes recurring. 

QtoD(x , flag := 0, a, b, c, d := 1111 , a^llst := []) := 

Prog 

If x = 0 

RETURN "O" 

[a := NUMERATOR (ABS (>0) , b := DENOMI NATOR (ABS (x ) ) ] 
c := APPEND(IF(x < 0, l, - n , l,n ), STRING (FLOOR (a, b))) 

Loop 

a := M0D(a, b) 

If a = 0 

If flag = 0 
If d = 1,11 
RETURN c 

RETURN APPEND(c, 11 . 11 , d) 

RETURN APPEND(a_l 1st , [a]) 

If MEMBER?(a, a_list) 

If flag = 0 

RETURN APPEND(c, 11 . 11 , INSERT (” ' 11 , d, P0SITI0N(a, a_l)st))) 
RETURN APPENDfXl-ist, [a]) 
a^list := APPEND(a_l ist , [a]) 

[a :* 10, d := APPEND(d, STRING(FL00R(a, b)))] 

In this case, it does not matter whether you Simplify or Approximate since Derive is returning a string. 
Try Simplifying QtoD(6858/5555) to obtain “1.2’3456”. In fact, in the Algebra window, the double 
quotes are not shown - so this looks like a number! 


Remarks 

(i) QtoD simplifies the recurring decimal format down to its shortest (simplest) form. For exam- 
ple, QtoD(DtoQ(“1 .201 2’01 2”)) Simplifies to “1 ,2’0120”. 

(ii) Setting flag = 1 in QtoD will display the sequence of remainders used in the calculation of the 
decimal expansion string. For example, Simplify QtoD(1/7,1 ) to display: 

[1, 3, 2, 6, 4, 5, 1], 

(iii) I never cease to be amazed at the calculating power of Derive. On my laptop Derive processed 
the 4293 symbols of QtoD(777/8777) in less than 2 seconds! 
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(iv) A useful instruction for investigations is: 


Dlistfn) := VECTOR QtoDfm ■ n) , m, 


Simplify D_list(1/7) - and note how the six recurring digits of 1/7 are repeated 
cyclically in the multiples of 1/7. 


3. Decimals to Quotients and Quotients to Decimals on the TI-89 


Having successfully composed and tested the above procedures in Derive 6 I thought of trying some- 
thing similar on my TI-89 Titanium calculator. Using 2 nd >VAR-LINK>f1 I first created a new folder 
DQ to store the three function procedures below. It is also useful to make DQ the default folder before 
entering or transferring the functions. The TI-89 does not appear to have a function to calculate the 
first position of a character in a string or member in a list. The function pos(a,b) is a DIY (Do It 
Yourself) version of this. 

:pos(a,b) igtocKx) 

: Func :Func 

: Local n,n :Local a,b,c,d,p,q,e 


pos(.a, 

Func 

Local n,n 
Dn: dintb^n 
While a^m i d ( b , n , 1 > 

If n>n Then: Return n 

Endlf 

n+l*n 

EndWhi le 

Return n 

EndFunc 


dtoqtxJ 

Func 

Local a,b,c,d,s 

dintx^+a 

If leftAx, 1 )=" Then 
_ l+s: a-l*a: right (x, a>*x 
Else: 1+s: Endlf 
dq\pos< n . M ; x)+Db 
dq^posC ,x)^c 

If c<a Then 

10 A (c-b)*a0 A (a-c)-l)+d 
expr ( r i ght ( x , a-c ) ) /d+d 
Else: :0^d: Endlf 
If c>2 Then 
expr < 1 eft < x , c- 1 )) -^c 
Else: Q+c: Endlf 
Return s*<c+d) 
EndFunc 


If x=G Then: Return "0" 
Endlf 



getHun(abs(x))+a 

? etDenon ( abs ( x ) ) -*b 
f x<0 Then: n -"+c: Endlf 

c&str i ng ( f 1 oor (. a/b ) ) 

Loop 

nod<a,b)-»a 
If a=Q Then 
If d=" 11 Then: Return c 
Else: Return c&V'&d 
Endlf: Endlf 
dqYpos< fa> , e)-*p 
If p<q Then 
left<d, q-p>+a 
righted, p}*b 

Return " . " &a& &b 

Endlf 

auqnent< fa> , e)+e 

q+l+q: 10*a+a 

dSfStr ing( floor (a^b) )-*d 

EndLoop: EndFunc 


Opting for functions allows the operations to be carried out on the HOME screen - this is useful for 
testing inverse properties of dtoq and qtod. Although the TI-89 is much slower and the display more 
limited you can still, for example, display the recurring decimal expansions of multiples of 1/7 by 
ENTERing seq(qtod(m/7),m,1 ,7). However, on the TI-89, there does not appear to be a method of 
setting up a flag variable with a default value - hence these functions only display the end results. 
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Recurring Expansions using General Number Bases in Derive 6 


In Derive (using Options>Mode Settings>...) it is possible to set the InputBase and OutputBase to 
any integer between 2 and 36, inclusive. To adapt the Derive 6 procedures (in sections 1 and 2 above) 
to work with these we first need a function to convert a base setting into a strict number output: 

C_B(b) := [IF Cb = Decimal, RETURN 10), IF(b = Binary, RETURN 2), IF(b = Octal, 


RETURN S), IF(b = Hexadecimal, RETURN 16), RETURN b] 

Then, in the local variable list of each procedure, add an assigned variable 

base := C_B(OutputBase) 

and within the procedure replace each occurrence of ”10” by ’’base”. 

This accommodates number bases < 10. 


For number bases > 10 (decimal) Derive uses letters A, B, C, ... for successive digits beyond 9 and, if 
the leading digit is >9, the number is prefixed by a zero. To follow this notation extend EXPR to: 

EXPR(x, n := 0, d, base := C_B(OutputBase)) := 

Prog 

x := NAME_T0_C0DES(x) 

Loop 

If x = [] 

RETURN n 

d := FIRST 00 - IF(FIRST(x) >64, 55, 48) 

If 0 < d < base 
n := base-n + d 
RETURN "DIGIT ERROR 1 ’ 
x := REST 00 


and alter the last line of QtoD to: 

[a :* base, d := APPEND(d, (STRING(FL00R(a, b)))i(-l))] 


Using alternative number bases can be confusing. In Derive, a straightforward method is to set the 
InputBase to Decimal and select the OutputBase. 

For example, with lnputBase:=Decimal and OutputBase:=Octal enter: DtoQ(QtoD(17/33)). 

f 21 T 

In the Algebra window this is displayed as: DtoQ QtoD 

l 41 Jj 

Highlighting and Simplifying QtoD... yields: DtoQ(0. '40760337017) . 

(This is the Octal recurring expansion of the Octal fraction!) 

21 

and Simplifying again gets back to: 

41 

In the reverse direction you have follow Derive notation when entering the string. 

For example, with OutputBase:=Hexadecimal enter: QtoD (DtoQ (“0 El. ABC ’123”)). 

In the Algebra window this is displayed as: Q to D ( DtoQ ( 0 El . ABC '123)) . 

r 4B348CCD 

Highlighting and Simplifying DtoQ. . . yields: QtoD 

l 555000 

(The fraction is Hexadecimal.) 

and Simplifying again gets back to: OE1 . ABC ' 123 . 
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Opportunities for experimentation are many. 

Try some activities using OutputBase:=Binary. 

(On the TI-89, to carry out similar activities using MODE settings of Binary or Hexadecimal, you 
would have to rewrite parts of the coding for the functions dtoq and qtod.) 

Working with these procedures I have extended my own experience of recurring decimals and recur- 
ring expansions well beyond hand-calculated examples. In my opinion, some form of precise recurring 
decimal/expansion notation ought to be a part of the number format of any IT product dedicated to 
mathematics and mathematical education. 

Reading Peter's contribution on periodic decimal numbers I remembered the procedure 
which I learned and applied in school time. See two examples: 


x = 3,53171 y = 0,23456 

lOOx = 353,171 171 lOy = 2,34563456 

lOOOOOx = 353171,171 171 lOOOOOy = 23456,3456 3456 

99900x = 352818 99990y = 23454 

x _ 352818 19601 _ 23454 1303 


99900 5550 99990 5555 

Now I tried to transfer this algorithm to DERIVE step by step. First I am doing it really step- 
wise and then I put all together to a function: 

#1: Notation := Decimal 

#2: PrecisionDigits := 20 

#3: NotationDigi ts := 20 

#4: [x := 3 . 53 '171, y := 0.2'3456, w := 10. '23] 

#5: codes(z) := NAME_T0_C0DES(z) 

Positions of decimal point and quotation mark 

#6: [dp(z) := POSITION (., z) , qp(z) := POSITION (', z)] 


dp(x) 

qp(x) 


' 2 

5 ‘ 

dp(y) 

qp(y) 

= 

2 

4 

. dp(w) 

qp(w) . 


. 3 

4 . 


The integer part of the number 

int_p(z) := C0DES_T0_NAME (VECTORC (codes (z) ) , i, dp(z) - 1)) 
#8: i 

#9: [i nt_p(x) , int_p(y), int_p(w)] = [3, 0, 10] 
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Which can be rewritten as 

per(z) 

fract(z) := tr_numb(z) + 

#24: pi (z) 

10 - 1 

' 19601 1303 1013 

#25: [fract(x), fract(y), fract(w)] = , , 

. 5550 5555 99 

I collect the whole procedure in a function: 

dn_to_f r(z , z_, dp, qp, ppl , pi, int_p, pp, per) := 

Prog 

[ z := NAME_T0_C0DES(z) , dp := POSITIONC" . " , z) , qp := POSITIONC" ' " , z)] 

i nt_p := CODES_TO_NAME(VECTOR(z_j,i , i , dp - 1)) 

#26: [ppl := qp - dp - 1, pi := DIM(z) - qp] 

pp := IFCppl = 0,0,CODES_TO_NAME(VECTOR(z_|i ,i ,dp + 1, qp - l))/10 A ppl) 
per := CODES_TO_NAME (VECTOR (z_ii , i ,qp + l,DIM(z)))/10 A (DIM(z) - dp - 1) 
int_p + pp + per + per/(10 A pl - 1) 

#27: [dn_to_f r(x) , dn_to_fr(y), dn_to_fr(w)] 

' 19601 1303 1013 

#28: , , 

. 5550 5555 99 

#29: [3.5317117117117117117, 0.23456345634563456345, 10.232323232323232323] 

Let's try a special case: 

#30: dn_to_f r(10 . '9) = 11 

I tested my function copying the decimal expansion of 777/8777 (= xx consisting of 4293 
numbers) from Peter's file and tried to convert it to a fraction 

7554973225475675059815426683377008089324370513842998746724393300672211461775093995670502449584140366- 

8679503247123162811894724849037256465762789107895636322205765067790816907827275834567619915688731912- 

9543124074285063233451065284265694428620257491170103680072917853480688162242223994531160988948387831- 

8332004101629258288709126125099692377805628346815540617523071664577873988834453685769625156659450837- 

4159735672781132505411871938019824541415062094109604648513159393870342941779651361513045459724279366- 

526147886521590520679047510538908510880710949071436709581861683946678819642246781360373703999 


777 

#32: dn_to_fr(xx) = 

8777 


It took 0.125 sec calculation time. (Peter's function DtoQ needs 0.234 sec). 

Is there also another way for the reverse task? 

Students might find out - analysing the calculation from page 22 - that the uncancelled de- 
nominator of the resulting fraction is always a number starting with a sequence of m nines 
followed by a sequence of n zeros. The number of zeros gives the length of the preperiod 
and the number of nines the length of the period. (Let them explain, why.) 

So we are looking for the minimum multiple of the given denominator which can be written as 
such a sequence of m nines and n zeros. 


My first task was to program an algorithm which delivers these possible denominators: 
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nines (k_, 1 := [9], 10, 11 := [9], a, aO, k, k_, i' , flag) := 

Pnog 

I := [ 9 ] 

II := [9] 
k := 1 
Loop 

10 := 11 

aO := FIRST(IO) -10 
11 := [aO] 

#33: Loop 

If DIM (10) = 0 exi t 
a := aO + FIRST(IO) 

11 := APPEND(11, [a]) 

10 := REST(10) 

1 := APPEND(1 , 11) 
k :+ 1 
If k = k_ 

RETURN 1 

#34: ni nes(8) = [9, 90, 99, 900, 990, 999, 9000, 9900, 9990, 9999, 90000, 
99000, 99900, 99990, 99999, 900000, 990000, 999000, 999900, 999990, 
999999, 9000000, 9900000, 9990000, 9999000, 9999900, 9999990, 9999999, 
90000000, 99000000, 99900000, 99990000, 99999000, 99999900, 99999990, 
99999999] 

Function fr_to_dn (see file Recurring_Josef, alas ... 

f 196Q1 ^ r 1303 'l r 1013 

#36: f r_to_d n , f r_to_dn , f r_to_dn 

L BBBO J L BBBB J L 99 

#37: [3 , B3 r 171, 0. 2 P 3456 , 10/23] 

f 93612037 

#38 : f r_to_dn = 1.1234B r 6789 

L S332BOOO J 

#39: P reci si onDi gits := 20 

#40: Notation Dig its := 20 

#41 : 1 . 12 34 B6 7896 7896 7896 7 

( 777 1 

#48 : f r_to_dn 

L S777 J 

#49 : 0 . r 088 3 268 3149139 79 715 7 2 2 00068 36 048 76 3814 BIB 2102 08499* 

SBB417B6S64B 3230033040902 3BS436S23B16007747B2193 232! 

B119 06UJS 26 364 36139911131366 07 04113 02 26 7 289 B 0666 B14 ' 

16 076108 009 B 7 0468 269 34 03 21294 291899 28 2 2148 79 7994 7 B9 1 
0462 B7263301S11BB 2922410846 B3070B2B2364133B30S19 186! 

B9 BS 7 B BS 3912498 B 7 BS 2 3174 2 0B 3 09 3 312 06 B6 26 06813 2619 34* 

B B 7 7 076449812 0086 BS99 B1008 31715 266 264 099 3 B 0B 7 B 36 74 3 ‘ 


. . . it works much slower than Peter's one (3900 sec for finding the correct result #49!) 


Some final comments: 

(1 ) It might be nice to check the derived formula #24 by manually calculation. 

(2) It is easy to extend the function for negative numbers. 

(3) I believe that this could be a fine problem for students - to analyse the procedure in 
general terms and then reproduce it using a CAS and finally write a function, which 
does the work. 


(4) Last but not least many thanks to Peter for this inspiring paper. Josef 
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Was verbirgt sich hinter Dr. Pest? 
What is hidden behind Dr Pest? 

Benno Grabinger, Neustadt/WeinstraBe, Germany 



Abbildung 1 


Der erste Griff am Morgen geht (oder sollte gehen) zur Zahnpastatube. Egal ob von Dr. Pest oder ei- 
nem anderen Hersteller, das Ding, das alleine auf seinem Verschluss stehen bleibt, verdient einige 
Beachtung. Industriell entsteht das Gebilde aus einem Zylinder, dessen ein Ende nach der Befiillung 
zusammengequetscht wird, so dass sich eine Falz ausbildet (Abbildung 2). 



Abbildung 2 


The first grip in the morning is (or should be) to the toothpaste tube. We would like to investi- 
gate the form of the tube and then calculate its volume. It is produced as a cylinder and its 
bottom is pressed together to a fold after filling it. Then we find the wellknown form. 


Can we compose this form using known geometric shapes? 
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Das Endprodukt liegt gut in der Hand, keine Ecken oder Kanten storen das glatte Erscheinungsbild. 


Lasst sich eine solche Tube aus bekannten geometrischen Grundkorpem zusammensetzen? 


Dazu muss eine Strecke AB, die Falz, mit den Punkten des Grundkreises der Tube verbunden werden. 
Da die breiteren Seiten der Tube sich eben anfuhlen, sucht man im Grundkreis diejenigen beiden 
Punkte P und Q, die eine zur Falz parallele Tangente besitzen. 


(Abbildung 3) 



The triangles APB and AQB are parts of the surface. The broad lateral faces of the tube feel 
nearly plane, so fold AB must be connected with those points of the base circle which have 
tangents parallel to AB. Connecting the points of the two semicircles with A and B we obtain 
to oblique circular cones. Figure 5 shows the two cones and figure 6 shows the complete 
figure including its cap. 


Die Dreiecke APB und AQB sind Teil der Tubenoberflache. 

Verbindet man nun noch jeden Punkt des Kreises mit A bzw. B, so entstehen zwei schiefe Halbkegel. 
(In der Abbildung 5 sind die zuvor gezeichneten Dreiecke wieder weggelassen, um die Kegel besser 
sehen zu konnen. Abbildung 6 zeigt die fertige Tube nebst Schraubverschluss.) 


Ich fuge die DERIVE-Konstruktion ein, da sie einen interessanten und wesentlichen Be- 
standteil dieses Beitrags darstellt. 

I include the DERIVE file, because the plotting procedure of the 3D-figure is an interesting 
and integrating part of this contribution. 

Josef 


D 


-N-L#70 


Benno Grabinger: What is hidden behind Dr Pest? 


p29 




Abbildung 5 


Abbildung 6 


Grundkreis und Strecke AB - Base Circle and segment AB 


#1: [1.5* SIN(t) , 1 . 5 -COS(t) , 0] 

#2: a := [-2.25, 0, 14.5] 

#3: b := [2.25, 0, 14.5] 

#4: [a, b] 



p und q sind die Kreispunkte in denen die Tangente parallel zur Strecke. 
p and q are points on the ci rcumference with tangents parallel to segment AB. 

#5: p := [0, 1.5, 0] 

#6: q := [0, -1.5, 0] 

Die vier Dreiecke im Raum werden erzeugt. 

The four triangles in space are generated. 
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#9: 


p q 

q a 


a p J 


# 10 : 


p q 
q b 


b p J 



k(t) sind die Punkte auf dem Grundkreis die einerseits mit a, andererseits mit b 
verbunden werden. 

Auf diese Weise entstehen 2 schiefe Halbkegel. 


k(t) are points of the base which are connected with a and with b. 
In this way two oblique cones are created. 

#11: k(t) := [1.5-COSCt), 1. 5 • SIN(t) , 0] 


#12 : VECTOR 


7T B 

[a, k(t) ] , t, — , tt , 0.1 

2 2 


#13 : VECTOR 


[b, k(t) ] , t, 


n 

2 


1 ' 

7T, 0.1 

2 


Der Verschluss - the cap 

#14: deckel := [1. 5-COS(t) , 1.5-SIN(t), s] 



End of the DERIVE file for plotting the figure. 


D 


-N-L#70 


Benno Grabinger: What is hidden behind Dr Pest? 


p31 


Welches Volumen besitzt die Tube? 

Das Volumen der beiden Halbkegel ist schnell angegeben. 1st r der Grundkreisradius und bezeichnet 

h die Hohe der Tube, so liefem die Halbkegel den Anteil 

, 11 , , 1 2 , 

2 r ■n-n = — -r -n-n 

2 3 3 


Der Restkorper besteht aus 4 Dreiecken, d.h. es handelt sich um eine auf einer ihrer Kanten stehenden 
Pyramide. (Abbildung 7) 


What is the volume of the tube? 

It is easy to find the volume of the two half cones. With 2 r = diameter of the base and h = 
height of the tube we receive 


2- — . — -r 2 -n-h 

2 3 



• n ■ h 


The remaining body is formed by four triangles and forms a pyramid standing on one of its 
edges (figure 7). 



Abbildung 7 Abbildung 8 

Die Berechnung desVolumens dieser Pyramide vereinfacht sich, wenn diese durch einen Schnitt in 
zwei gleich groBe Teile zerlegt wird. Die Schnittebene geht durch die Punkte A, B und M. In der Ab- 
bildung 8 wurde das vordere begrenzende Dreieck weggelassen um die Unterteilung der Pyramide in 
die zwei Teile zu sehen. 


For calculating the volume of this pyramid we divide it into two equal parts by an intersection 
plane ABM (which can be seen gray coloured in figure 8). ABM is base of pyramid PMBA. 

1 

Its height is half diameter of the tube r = PM. Its area is — f bwith f= AB. 


The volume of this solid is given by 2 • 


■fh 


j 


r = — f ■ h ■ r and the total volume 
3 


V =--r 2 -n-h + --f -h-r. 

3 3 
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Assuming that we fold the tube so that f equals half of the perimeter of the cylinder = r- i.e. 
1 

f = — -2- r ■ n = r -n - we obtain the final value for V\ 

2 


1 2 . 1 . 2 2 
V= — -r-n-h + — -r-n-h-r = — -r-n-h. 


1 

3 3 3 

The volume of the tube is 2/3 of the volume of the starting cylinder. 


The tube shown in figure 1 is given by r = 1.5 cm and h = 14.5 cm. Hence it is filled with 
2 

— -1,5 2 -Ti -14,5 « 68 cm 3 (ml) toothpaste. The label promises 75 ml. A similar difference can 
3 

be observed at other tubes. It might be explained by a bump out of the tube (which can also 
be found at milk boxes). The missing ml are hidden in these bumps. 


Damit ist das Volumen des Restkorpers gleich 2 • 


• f'h \-r — — ■ f -h-r . 


1 2 1 

Fur das Tubenvolumen V ergibt sich dann zu: V - — -r -n-h-\ f -h-r 

3 3 

Nimmt man an, dass der Falzvorgang so stattfindet, dass der halbe Zylinderumfang gleich / wird, d.h. 

/ = --2-r-n = r-n, dann gilt fur V: 

2 


T r 1 2 7 1 7 2 2 ; 

V — — -r - Ti • h -\ r -n- h-r = — -r -n-h 


. 2 


Das Tubenvolumen betragt damit — des Volumens des urspriinglichen Zylinders der befiillt wird. 


Fur die Tube aus Abbildung 1 ergeben sich beim Abmessen die folgenden Werte: 
r = 1 ,5 cm und h= 14,5 cm. 

2 2 3 

Die Rechnung liefert dann: — - 1,5 • n - 14, 5 « 68 cm (ml) Zahnpasta. 

Die Tubenaufschrift verspricht dagegen 75 ml Inhalt. Dieser Effekt ist auch bei weiteren untersuchten 
Tuben zu beobachten. Eine mogliche Erklarung ist eine Ausbeulung der Tube wie dies auch von 
Milchtiiten bekannt ist. In dieser Ausbeulung versteckt sich dann das fehlende Volumen. 


Samtliche Abbildungen wurden mit dem Programm DERIVE angefertigt. Fur Hinweise zur Volumen- 
berechnung bedanke ich mich bei Herm Dr. Chr. Fahse aus Neustadt. Fur die technische Unterstiit- 
zung bedanke ich mich bei Dr. Klaus Wagner aus Neustadt. 

All figures were produced by using DERIVE. I am very grateful for advice to the calculation of 
te volume (Dr Chr Fahse, Neustadt) and for Technical support (Dr. Klaus Wagner, Neustadt). 
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I made a typo in surface #5 from DNL#69 and received #5a: 

Surface #5a: {^y 2 +y 2 j =x 2 y 2 (^z 2 + lj 



Autograph MUPAD ] 


[1] Produced using i mpl i ci t_peter . dfw from DNL#64 

6 2 2 2 2 2 2 

ImplicitDots(8*y = x -y • (z + 1) a x +y +z < 25, [-5, -5, -5], [5, 5, 5], 0.2). 

[2] You can download Surfer for free at www . imaginary2008 . de. The surface looks a bit strange 
because the Surfer plots are bounded by a sphere. 

Surface #6: ^ 2 -z 3 j =^x+ y 2 ^- z 3 - The UFO 



Surfer 


DPGraph 
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derive [3] 


Surface #7: y 2 +z 3 =z 4 + x 2 z 2 


Autograph 



DERIVE 


DPGraph 


2 3 4 2 2 

VECTOR(ContourPts_XY(y + z - z + x *z , 


- 3 , 3 , - 3 , 3 , 0 . 2 , 0 . 2 ), 1 , - 3 , 3 , 0 . 2 ) 



Autograph Surfer 

[3] The next DERIVE plots are produced by using the plot routines from pol ycontou r . dfw (presented 
in DNL#63). They are included in file FAZ5 . dfw 
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pip2 + plql + plq2 - pi 
plp2 

- + p2ql + p2q2 - p2 

Rplp2 

plql p2ql 

-i- - qlq 2 + ql 

#4: Rplql Rp2ql 

p2q2 plq2 qlq2 

+ + + q2 

Rp2q2 Rplq2 Rqlq2 

aB-plp2 + bB-plql + eB-p2q2 

f6-qlq2 + c6-plq2 + d6-p2ql 

DERIVE does not disclose the solution. Artur sent another equation with only four variables - 
together with the MATHEMATICA-solution. 

#1: CaseMode := Sensi tive 

#2: InputMode := Word 

Solve for plp2, plq2, p2ql and qlq2. 

I" Rplql- (Rp2ql- 


Rplp2- (Rplq2- (Rp2q2- (Rqlq2-q2 + ■ 


#3: 

LAI- (LA2 ■ Rplql- (Rp2q2 - 1). (Rplp2- (p2 - p2ql) + plp2) + Rp2q2- (Rplp2- (Rplql- (pi - plq2) - pi + plp2 


LB1- (LB2 ■ Rplq2 ■ Rqlq2 ■ p2ql- (Rp2ql - 1) + Rp2ql- (Rplq2 ■ (Rqlq2 ■ (plq2 + qlq 2) - ql 


Here are the MATHEMATICA solutions pi p2_, pi q2_, p2q1_and q1q2 


#4 : plp2_ = 

LAI- (LA2- (Rp2q2 - 1) ■ (LB1- (LB2- (Rp2ql - 1) ■ (Rplql- Rqlq2 - Rplq2) + (Rp0q2 - Rqlq2) ■ (Rplp2 ■ Rp2q! 

Rplq2- (LAI- (LA2- (Rp2q2 - 1) ■ (LB1- (LB2- (Rp2ql - l)-(Rplp2-p2 + Rplql- (Rqlq2-q2 + ql) + pi) + . 

#B : plq2_ = 

LAI- (LA2- (Rp2q2 - 1) ■ (LB1- (LB2- (Rp2ql - 1) ■ (Rplql- Rqlq2 - Rplq2) 


(1 ) I remove the denominators and obtain four linear equations: 

#8: eql := Rplql ■ (Rp2ql ■ (ql - qlq2) + p2ql) + Rp2ql'(pl - plp2 - plq2) 

#9: eq2 := Rplp2 ■ (Rplq2 ■ (Rp2q2 ■ (Rqlq2 ■ q2 + qlq2) + Rqlq2'(p2 - p2ql)) + Rp2q2 ■ Rqlq2 ■ p 

#10: eqi LAI ■ (LA2 ■ Rplql ■ (Rp2q2 - 1) ■ (Rplp2 ■ (p2 - p2ql) + plp2) + Rp2q2 ■ (Rplp2 ■ (Rplc 
LA2 • Rplql ■ (plp2 - Rplp2 ■ f Rp2q2 • f plp2 + p2 - p2ql) - p2 + p2ql)) + Rplp2 ■ Rp2q2 ■ 


#11: eq4 := LB1 ■ (LB2 • Rplq2 ■ Rqlq2 ■ p2ql • (Rp2ql - 1) + Rp2ql ■ CRplq2 • (Rqlq2 • (plq2 + qlq2) 

Rp2ql ■ (Rqlq2 ■ (p2ql + qlq2) - qlq2)) - Rp2ql ■ Rqlq2 ■ plq2 ■ (Rplq2 - 1) 

The expressions are so large that I show only part of the screen. The file is among 
mth70.zip. 

DERIVE is also unable to solve this equation applying SOLVE or SOLUTIONS. I tried to fool 
my old friend DERIVE. 
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(2) In order to have an easier reading I substitute for pi p2, pi q2, p2q1 and ql q2 the varial 

#12: eql := SUBST(eql, [p0_p2 , pOr^, p2ql, qlq2], [w, x, y, z]) 

#13: eql := - Rp2ql-x + Rplql- y - Rp2ql- (Rplql- z - Rplql- ql - pi + w) 

#14: eq2 := SUBST(eq2, [plp2, plq2 l p2ql, qlq2], [w, x, y, z]) 

#15: eq2 := Rplp2- Rp2q2- Rqlq2-x - Rplq2- (Rplp2- Rqlq2-y - Rplp2- Rp2q2- z - Rplp2- Rqlq2- ( 

#16: eq3 := SUBST(eq3, [plp2, plq2, p2ql, qlq2], [w, x, y, z]) 

#17: eq3 := - x ■ (LAI- Rplp2 ■ Rp2q2 ■ (Rplql - 1 ) + Rplp2- Rp2q2- (1 - Rplql)) - y-(LAl-LA2-R 

+ LAI- (LA2- Rplql- (Rp2q2 - l)-(Rplp2-p2 + w) + Rp2q2- (Rplp2- (Rplql- pi - pi + w) 
p2)) + Rplp2- Rp2q2- (Rplql - 1) ■ (w - pi) 

#18: eq4 := SUBST(eq4, [plp2, plq2, p2ql, qlq2], [w, x, y, z]) 

#19 : eq4 := x- (LB1- Rp2ql- Rqlq2- (Rplq2 - 1) - Rp2ql- Rqlq2- (Rplq2 - 1)) + Rplq2- (LB2- Rql 

LB2) ■ (Rqlq2 - 1)) 

(3) The system is of the form: a11w + a12x + a13y + a14z + a10 = 0 

a21 w + a22 x + a23 y + a24 z + a20 = 0 

a31 w + a32 x + a33 y + a34 z + a30 = 0 

a41 w + a42 x + a43 y + a44 z + a40 = 0 

and I determine all coefficients, starting with al 1 , a21 , a31 , a41 in one step: 

d 

#20: coll := — [eql, eq2, eq3, eq4] 

dw 

d 

#21: col 2 := — [eql, eq2 , eq3 , eq4] 

dx 

d 

#22: co!3 := — [eql, eq2 , eq3 , eq4] 

dy 

d 

#23: co!4 := — [eql, eq2 , eq3 , eq4] 

dz 

#24: const := SUBST([eql, eq2 , eq3 , eq4], [w, x, y, z], [0, 0, 0, 0]) 

I try applying ROW_REDUCE 

#25: sys_mat := [coll, co!2, co!3, co!4, -const]’ 

#26 : RQW_R EDUCE (sy s_mat) 

I cancel calculation - no result after waiting a long time. 

Another try: 

sols := (5QLUTIQNS([aU_- w + al2-x + al3-y + al4-z + alQ = 0, a21-w + a22-x + a23»- 

# 27 : 

a30 = 0, a41-w + a42-x ■+■ a43-y + a44-z + a40 = 0], [w, x, y, z])) 

1 
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plp2 := sols , plq2 := sols , p2ql := sols , qlq2 := sols 1 
#2S : L 1 2 3 4 J 


This is the general solution. Now I resubstitute for the coefficients: 


#29: 

[all 

:= coll , 
1 

a21 

:= coll , 
2 

a31 :z coll , 
3 

a41 coll , 
4 

alO := const 

1 

#30: 

[al2 

:= co!2 , 
1 

a22 

:= co!2 , 
2 

a32 := co!2 , 
3 

a42 := co!2 , 
4 

a20 := const 

2 

#31: 

[al3 

:= co!3 , 
1 

a23 

:= col 3 , 
2 

a33 := co!3 , 
3 

a43 := co!3 , 
4 

a30 := const 

3 

#32: 

[al4 

co!4 , 
1 

a24 

co!4 , 
2 

a34 :z co!4 , 
3 

a44 co!4 , 
4 

a40 const 

4 


This should be the solution: 

#34: plp2 = 


LAl-(LA2-(Rp2q2 - 1) ■ (LB1- (LB2- (Rp2ql - 1) ■ (Rplql- Rqlq2 - Rplq2) + (Rp0jq2 - Rq 


- 1) + Rplql- (1 - Rp2ql-Rqlq2) + Rplq2-(Rp2ql - 1)) + (1 - Rplq2) ■ (Rplp2- Rp2ql 


- Rqlq2) + Rplql- Rqlq2- (Rp2q2 - Rplq2) + Rplq2- (Rqlq2 - Rp2q2)) + Rplp2-(Rplq2 
Rp2ql - 1)- (Rplql- Rqlq2 - Rplq2) + Rplql- (Rqlq2- (Rp2ql- (ql + q2) + p2) - Rp0jq2- 

- Rp2q2 - Rqlq2 + 1) - Rqlq2- (Rp2ql - Rp2q2)) + Rplq2- (1 - Rqlq2) ■ (Rp2ql - 1)) 
2 - pi - q2) + Rqlq2-pl)) - LB2- (Rplql- (Rp2ql- (Rqlq2- (p2 + ql) - ql) - p2) + Rp 
q2) - Rp2q2- (Rp2ql-Rqlq2 - 1)) + Rplql-Cl - Rp2ql- Rqlq2) ■ (Rplq2 - Rp2q2) + Rplq 

and so on for all unkowns ... and finally: 

I compare your_solutions from #4 to #7 with my solutions: 

#37: [qlq2 - qlq2_, pip2 - plp2_, p2ql - p2ql_, plq2 - plq2_] 

#3S : [0, 0, 0, 0] 

The solutions are identical. 

So this procedure could be a model howto tackle the more complicated system from the other file. 


I treid to solve Artur's equations using the open source CAS Maxima and succeeded with my 
first attempts. Calculation time was incredible short. See parts of the output on the next 
page. 

We - Artur and I - are wondering if anybody will be able to solve the given simultaneous 
equations in a satisfying way using DERIVE. 
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Here are screen shots from the wxMaxima - output. 

See first the example with four unknowns: 

£%i4) eq4 : LBI* (LB2 *Rplq2 *Rqlq2 *p2ql* (Rp2ql - 1) + Rp2ql* (Rplq2 * (Rqlq2 * (plq2 + qlq2 ) - qlq2 ) 

- Rqlq2*plq2)) + LB2 *Rp lq2 * ( Rqlq2 *p2 ql - Rp2ql* (Rqlq2 * (p2ql + qlq2 ) - qlq2 ) ) - Rp2ql*R( 

- i) ; 

(%o4) LBI ( p2ql Rplq2 [ Rp2ql - 1 ) Rqlq2 LB2 + Rp2ql [ Rplq2 [ ( qlq2 + plq2) Rqlq2 - qlq2) - plq2 Rqlq2 ) ) + Rp. 
[ p2ql Rqlq2 - Rp2ql (( qlq2 + p2ql ) Rqlq2 - qlq2 )) L B2 - plq2[ Rplq2 - 1 ) Rp2ql Rqlq2 


(%i5) sol : solve ( [eql, eq2 , eq3 , eq4] , [plp2 , plq2, p2ql, qlq2] ) ; 

« Ausciruck zu lanq u m anqezeiqt zu werden! » 

(% 16) sol [1] [1] ; 

(%o6) plp2 = - ( Rp2ql ( Rplql { Rqlq2 ( Rp2q2 ( Rplp2{ p2 LA2 ( ( LEI - 1 ) LB2 + LAI ( 1 - LB1 ) LB2 ) + q2 LA2 ( Lj 
Rplp2 Rplq2 { ( LEI - 1 ) LE2 + LAI { 1 - LEI ) LE2 + LA2 ( - LEI + LAI ( LEI - 1 ) + 1 ) ) ; + Rplp2 
( p2 LA2 ( LAI ( LB1 - 1 ) LE2 + ( 1 - LEI ) LE2 ) + q2 LA2 { - LB1 + LAI [ LB1 - 1 ) + 1 ) ) + Rplp2 Rpl q2 
{ p2 { ( LEI - 1 ) LB2 + LAI { 1 - LSI ) LB2 ) + q2 LA2 ( LEI + LAI { 1 - LSI ) - 1 ) ) ; + ql ( Rql q2 

{ Rplp2 Rp2q2 LA2 { LAI { LB2 - LB1 ) - LB2 + LEI ) + Rplp2 Rplq2 { LAI { LB2 - LB1 ) - LB2 + LB1 ) + Rplp2 LA2 { LB2 + L 
( Rplp2 Rplq2 ( ( LSI - 1 ) LB2 + LAI ( 1 - LEI ) LB2 + LA2 { - LEI + LAI ( LSI - 1 ) + 1 ) ) + Rplp2 LA2 ( LB2 + LAI ( 1 - LBi 
[ LB2 + LAI [ LBI - LB2 ) + LA2 [ LBI + LAI [ 1 - LBI ) - 1 ) - LBI ) + Rplp2 LA2 [ - LB2 + LAI [ LB2 - 1 ) + 1 ))+ pi Rqlq2 


This is the second one - using the form with generalized coefficients (expression #4 from the 
above DERIVE file): 

(% i5) eqS : a5*plp2+b5*plql+e5*p2q2=0; 

(%o5) e5 p2q2 + bS plql + a5 plp2 = 0 

(% i 6 ) eq6 : f 6*qlq2+c6*plq2+d6*p2ql = 0; 

(%o6) f6qlq2+d6 p2ql + c6 plq2 = 0 

(% i7) sol : solve ([eql,eq2,eq3,eq i l,eq5,eq6] , [plp2,plql,plq2,p2ql,p2q2,qlq2] ) ; 

(%o7) [ [ plp2 = { Rplp2 ( Rp2ql (Rplql 

( ql ( ( e5 f€ Rp2q2 + b5 f6 Rplq2 ) Rqlq2 +[ b5 d6 - c6 e5 ) Rplq2 Rp2q2 )+[ b5 d6 - c6 e5 ) q2 Rplq2 Rp2q2 
p2 ( Rplql (( e5 f6 Rp2q2 + b5 f€ Rplq2) Rqlq2 - c6 e5 Rplq2 Rp2q2) + Rp2ql ( Rplql ( e5 Rp2q2 + b5 d6 
f€ q2 Rplql Rplq2 Rp2q2 Rqlq2 ) + pi Rplp2 

( Rp2ql ( Rplql ( b5 d6 Rp2q2 + b& c6 Rplq2) Rqlq2 + e5 f€ Rp2q2 Rqlq2 - c6 e& Rplq2 Rp2q2) + Rplql ( b5 
( Rplp2 ( Rp2ql ( ( e5 f6 Rp2q2 + a5 f€ Rplq2) Rqlq2 + Rplql [[ b5 d6 - a5 d6) Rp2q2 + ( bS c6 - a5 c6) Rpl 
+ Rplql ( ( b5 f€ - a5 f€) Rp2q2 Rqlq2 + ( aS c€ - b5 c€) Rplq2 Rp2q2 ) ) + Rplql ( ( - e5 f€ Rp2q2 - b5 f€ R, 

( Rplql ( - d6 e5 Rp2q2 - b5 d6 Rplq2) Rqlq2 + d6 e5 Rplq2 Rp2q2) ) , plql = - ( Rplp2 ( Rplql Rp2ql 
[ ql ( ( e5 f6 Rp2q2 + a5 f6 Rplq2 ) Rqlq2 +[ a5 d6 - c6 e5 ) Rplq2 Rp2q2 )+[ a5 d6 - c6 e5 ) q2 Rplq2 Rp2q2 
( Splql ( ( e5 f€ Rp2q2 + a5 ft Rplq2) Rqlq2 - c6 e5 Rplq2 Rp2q2) + Rplql Rp2ql ( e5 Rp2q2 + a5 Rpl 
Rqlq2 )+ pi ( Rplp2[ Rplql ( a5 f€ Rp2q2 Rqlq2 - a5 c6 Rplq2 Rp2q2) + Rplql Rp2ql ( a5 Rp2q2 + a5 c 
( e5 f€ Rp2q2 Rqlq2 - c6 e5 Rplq2 Rp2q2) + ci^ e5 Rplql Rp2ql Rp2q2 Rqlq2 ) + Rplql Rp2ql ( e5 q2 Rp 
+ e5 f& q2 Rplql Rplq2 Rp2q2 Rqlq2 ) / ( Rplp2 ( Rp2ql 
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Titbits from Algebra and Number Theory (35) 
or Yet Another Treatise on RSA 

(c) Johann Wiesenbauer, Vienna University of Technology 


Hm, RSA again - not exactly very imaginative you might think. In fact, a huge amount of 
articles are referring to it, as a simple Google search proves. Isn't everything said and 
done, when it comes to RSA? Well, I will come back to this question before long. Let me 
point out first that I'm focussing here mostly on the mathematical background of RSA 
and the implementation of certain algorithms rather than the cryptographical environ- 
ment. In particular, I will assume here a basic knowledge what RSA is all about (if 
needed, an excellent reference for cryptography is 

http://www.cacr.math.uwaterloo.ca/hac/') . 

Basically, if n is a natural number that is a product of two "big" primes p and q, and x is 
a message coded as an element of the residue class ring Zn, i.e. x £ {0,l,2,..,n-l}, then 
the encryption E is a bijective mapping from Zn onto Zn of the form x -> x A e mod n for 
some fixed natural number e. Now it's a simple algebraic exercise to prove that two 
exponents el and e2 induce the same mapping if and only if 

el=e2 mod A(n), where A(n)=lcm(p-l,q-l). 

In particular, we may assume w.I.o.g. that e < A(n). This innocuous looking assertion has 
the important implication that the mapping induced by the exponent el is the identity 
(represented by the exponent e2=l), if and only if el = 1 mod A(n), which in turn implies 
that the mappings belonging to the exponents e and d are inverse to each other if and 
only if 

de = 1 mod A(n). 

For some of you this might be the first surprise as you might have been used to the 
condition 

de = 1 mod q>(n), where q>(n)=(p-l)(q-l). 

Furthermore, in the "classical" RSA the exponents e and d were supposed to be < q>(n) 
rather than < A(n). Well, as we know, it works too, because A(n) is a divisor of (p(n), but 
the exponent d is usually not best possible as to its size. Let's do a small example, just 
to see what can happen. We use n= 15251=101*151 for this, and compute d for e=301, 
using the classical RSA-scheme. Since A(n)= lcm(100,150)=300 and 301=1 mod 300, we 
know that the mapping belonging to exponent e is the identity mapping. Hence, we would 
d expect to be 1. Unfortunately this not true. 
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#1: d := INVERSE_MOD(301, EULER_PHI(15251)) 

#2: d := 14701 


The number d looks quite inconspicuous, doesn't it? Nevertheless, the corresponding 
mapping is actually the identity!!! 

14701 

#3: EVERY(MOD(x , 15251) = x, x, 0, 15250) = true 

Admittedly, this is the worst case, and usually d is only a few bits larger than neces- 
sary, but all the same is clearly better, as it is optimal, when it comes to size, and its 
computation is still very simple. Sadly enough, even though, still rather few authors use 
X(n), one of them is 

http://www.staff.uni-mainz.de/pommeren/Kryptoloqie/Asvmmetrisch/ 

(Sorry, as it is only for people with a basic command of German!) 


This suggests the following question: What is general relationship between <p(n) and A(n) 
for any positive integer n? Well, <p(n) is the number of elements of Zn* = (a | l<=a<=n and 
gcd(a,n)=l}, and <p is called Euler's ^-function, after its first investigator, whereas A(n) 
is the smallest positive integer k such that a / 'k= 1 mod n for all a in Zn*. In short, q>(n) is 
the order and A(n) the so-called exponent of Zn*. <p(n) is the library function 
euler_phi(n), which we already used above, and A(n) , which is usually called after Carmi- 
chael, can be implemented as follows. 

carmichael_A(n) := 

Prog 

#4: If M0D(n , 8) = 0 

n :/ 2 

LCM(VECTOR(f_4,l A f_j,2 • (1 - l/f_|l) , f_, FACTORS(n))) 

#5: carmichael_A(15251) = 300 

This raises the following interesting question: Does our RSA-scheme above work for an 
arbitrary n using this definition of A(n) ? Unfortunately, the ans\Nzr is no in general, 
though it does work if n is squarefree, i.e. the product of different prime numbers, 
which is at least a slight generalization of the original RSA, where n is the product of 2 
different primes. On the other hand, even if n is not squarefree, it still works for 
"most" x, or more precisely it works, unless gcd(x,n/gcd(x,n)) > 1. 


Let's test this for n=45, hence A(n)=12, and e=d=5. 
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5 5 

#6: SELECT(MOD((x ) , 45) * x, x, 0, 44) = [3, 6, 12, 15, 21, 24, 30, 33, 39, 42] 

45 


#7: SELECT|GCD[ 

v v GCD(x , 45) 


( 45 

) x, 

l GCD(x , 45) 


> 1, x, 0, 44 


[3, 6, 12, 15, 21, 24, 30, 33, 


39, 42] 


However, as for all these "generalizations" of RSA, we are not so much concerned 
about those exceptional values of x - for an n with "big" prime factors the probability 
to encounter any of those is virtually zero but about the fact that the factorization 
problem becomes far easier, if n splits up into many primes! And anyone, who can factor 
n, can compute A(n) easily (see the algorithm #4 above!) and a fortiori the secret key d, 
which is the inverse of e mod A(n). 

Let's turn to a realistic RSA-implementation now. To this end I will provide a routine 
setup(b,e) that sets a number of global variables, namely 

- the primes p,q with p t q of about the size 2^(b/2) 

- their product n = p*q, whose binary representation has exactly b bits, where b=1024 
by default 

- the public exponent e, where e =2*16+1 by default 

- the corresponding private exponent d = e*(-l) mod A(n) 

- the values dp and dq, which are d mod (p-1) and d mod (q-1) by definition, respectively, 
and will be used for decryption later on 

setup(b := 1024, eO := 2 A 16 + 1) := 

Prog 
e := eO 
Loop 

p := NEXT_PRIME(RANDOM(2 A FLOOR(b, 2))) 

If GCD(e, p - 1) = 1 exit 
Loop 

#8: n := 2 A (b - 1) + RAND0M(2 A (b - 1)) 

q := NEXT_PRIME(FLOOR(n , p)) 
n := p-q 

d := INVERSE_MOD(e , LCM(p - 1, q-1)) 
dp := MOD(d, p-1) 
dq := MOD(d, q-1) 

If p t q a FLOOR(LOG(n , 2)) = b - 1 a NUMBER?(d) exit 

The following call of setup() that provides all global variables above for the default val- 
ues of b=1024 and e=2*16+l takes usually only tenths of a second. 


#9: 


setupQ = true 
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Now, let's assume that the message m is given as string of symbols from the extended 
ASCII-character set. In a first step, we concatenate all 8-bit representations of the 
characters of m, getting thus a single binary string s. Next we split up an expanded ver- 
sion of s into k binary blocks of the fixed length II for some II < b. (We will choose the 
maximal possible value ll=b-l in our implementation below.) In order to get that ex- 
panded version of s mentioned above, we first adjoin the bit 1 at the rightmost position, 
and thereafter the smallest number of zeros such that the length of the new s becomes 
divisible by I (zeropadding). Now after deleting leading zeros these k binary blocks of 
length I are interpreted as binary representations of nonnegative integers xl,x2,...,xk all 
below n, which are encrypted then into yl,y2,..,y2, using the RSA-scheme at issue. These 
y-values are usually converted into binary blocks of some fixed length 12, where I2>=b, 
and concatenated. For didactic reasons I decided to use the explicit vectors 
(xl,x2,..,xk) and (yl,x2,...,yk) of decimal values in the following though. In the following, 
I wrote separate routines for the source encoding/decoding and the actual encrypt- 
ing/decrypting, respectively. Furthermore I give two versions of decrypt(c), where the 
faster second one makes use of the Chinese Remainder Theorem (CRT) and the global 
variables dp and dq, we provided in the setup above. 

encoding(t, b_, c_ := 1, s_ 0) 

Prog 

b_ := FLOOR(LOG(n , 2)) 
t := NAME_TO_CODES(t) 

Loop 

If t = [] exit 
c_ :+ 8 

s_ SHIFT(s_, 8) + FIRST(t) 
t := REST(t) 

#10: s_ := 2-s_ + 1 

s_ := SHIFT(s_, MOD(-c_, b_)) 
t := [] 

c_ SHIFTCl, b_) 

Loop 

If s_ = 0 
RETURN t 

t := ADIOIN(MOD(s_, c_), t) 
s_ := SHIFT(s_, -b_) 


#11: encrypt(t) := MAP_LIST(MOD(x , n) , x, encoding(t)) 

decoding(v, b_, s_ := 0, t_ := "") := 

Prog 

b_ := FLOOR(LOG(n , 2)) 

Loop 

If v = [] exit 

s_ := SHIFT(s_, b_) + FIRST(v) 
v := REST(v) 

Loop 

#12: b_ := MOD(s_, 2) 

s_ := SHIFT(s_, -1) 

If b_ = 1 exit 
Loop 

If s_ = 0 
RETURN t_ 

t_ := AD30IN(C0DES_T0_NAME(M0D(s_, 256)), t_) 
s_ := SHIFT(s_, -8) 
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d 

#13: decrypt(c) := decodi ng(MAP_LIST(MOD(y , n), y, c)) 

[dp dq 1 

#14: decrypt(c) := decoding(MAP_LIST(CRT(LMOD(y , p) , MOD(y , q)J, [p, q]), y, c)) 

Let's check it using a famous piece of text from a speech of Abraham Lincoln ("Gettys- 
burg address") consisting of 1494 characters. I also used here for the very first time a 
timing function that will give you the run time in seconds on the screen itself rather 
than on the right side of the status bar. (See the examples below as to its use.) In par- 
ticular you can see that the total time for both encrypting and decrypting was only 
0.48s on my PC and most of this time, namely 0.41s, was spent on decrypting a vector 
with 12 numbers each with about 308-309 digits. Hence, the RSA decryption took about 
0.034s per block. 

timing(t, r) := 

Prog 

t APPEND (STRING ( RANDOM (0) - t) , "s") 

#15: Loop 

If DIM(t) < 4 

t := ADJOIN("0" , t) 

RETURN INSERT(" . " , t, -3) 

#16: t := November 19, 1863 — Four score and seven years ago, our fathers brought 

forth upon this continent a new nation: conceived in liberty, and dedicated 
to the proposition that all men are created equal. Now we are engaged in a 
great civil war . . . testi ng whether that nation, or any nation so conceived 
and so dedicated... can long endure. We are met on a great battlefield of 
that war. We have come to dedicate a portion of that field as a final 
resting place for those who here gave their lives that that nation might 
live. It is altogether fitting and proper that we should do this. But, in a 
larger sense, we cannot dedi cate . . . we cannot consecrate... we cannot hallow 
this ground. The brave men, living and dead, who struggled here have 
consecrated it, far above our poor power to add or detract. The world will 
little note, nor long remember, what we say here, but it can never forget 
what they did here. It is for us the living, rather, to be dedicated here 
to the unfinished work which they who fought here have thus far so nobly 
advanced. It is rather for us to be here dedicated to the great task 
remaining before us... that from these honored dead we take increased 
devotion to that cause for which they gave the last full measure of 
devotion... that we here highly resolve that these dead shall not have died 
in vain... that this nation, under God, shall have a new birth of freedom... 
and that government of the people... by the people... for the people... shall 
not perish from the earth. 

#17: [DIM(t) , DIM(n) ] = [1494, 309] 

#18: [DIM(encodi ng(t)) , timi ng(RANDOM(0) , encodi ng(t)) ] = [12, 0.07s] 

#19: [DIM(encrypt(t)) , timi ng(RANDOM(0) , encrypt(t))] = [12, 0.08s] 

#20: DIM(n) = 309 

#21: S0LVE(decrypt(encrypt(t)) = t) = true 

timi ng(RANDOM(0) , decrypt(encrypt(t))) = 0.48s 


# 22 : 
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Now let's discuss some security aspects of RSA. In the first place, the primality tests 
to find the primes p and q in our routine setup() are only probabilistic, i.e. there is a tiny 
chance that one of them is actually composite. This raises the interesting question: Is it 
possible that our RSA-encryption and decryption routines above still work, if one of the 
numbers p and q is not prime or even both? Surprisingly enough, the answer is yes. 

To be more precise, there are actually 3 properties we really need regarding the posi- 
tive integers p and q: 

1) a*p = a mod p for all integers a 

2) a~p = a mod q for all integers a 

3) gcd(p,q) = 1 

An integer n>l that fulfils a*n=a for all integers a without being a prime is called a 
Carmichael number. Unfortunately, those numbers exist and there are even infinitely 
many of them (http://www.math.dartmouth.edu/~carlp/PbF/paper95.pdf ) Due to a 
theorem by Korselt, they can be characterized as squarefree integers n > 1 such p - 1 
divides n - 1 for any prime divisor p of n. It can be easily proven then that a Carmichael 
number must always be odd and has at least 3 different prime factors. Furthermore, 
any n of the form n=(6k+l)(12k+l)(18k+l) for some positive integer k, where 6k+l,12k+l 
and 18k+l are prime, fullfills the conditions of Korselt's theorem. 

We use the latter fact to construct Carmichael numbers greater or equal to a given 
number s, which are used then in an RSA-emulation for out text t above that works in- 
deed. Needless to say that this is a horror scenario in practice, as the factoring prob- 
lem for n and hence the computation of A(n) and d becomes so much easier for an at- 
tacker. Note also the relatively small value of A(n), which forces d to be rather small, 
too! 


Carmichael ?(n) := 

Prog 

If PRIME?(n) v (n - l).MOD(n, 2) = 0 
#23: RETURN false 

If SOME(e_ > 1, e_, (FACTORS(n)) COL 2) 

RETURN false 

EVERY (MOD (n - 1, p_ - 1) = 0, p_, ( FACTORS (n)) COL 1) 

#24: ITERATE(IF(PRIME(6-k + 1) a PRIME(12-k + 1) a PRIME(18-k +1), k, k + 1) , 

k, 1000) = 1025 
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#25: NS0LVE((6-k + l).(12-k + l)-(18-k + 1) = 1000, k) = (k = -0.56104944 - 

0.7932487967 - L v k = -0.56104944 + 0.7932487967 - L v k = 0.8165433244) 

#26: NSOLVE((6- k + l)-(12-k + l)-(18-k + 1) = 1000, k, Real) = (k = 0.8165433244) 

find_CN(s, k_) := 

Prog 

k_ := CEILING(RHS(NS0LVE((6- k + l)-(12-k + l)-(18-k + 1) = s, k, Real))) 

#27: Loop 

If PRIME(6-k_ + 1) a PRIME(12 • k_ + 1) a PRIME(18-k_ + 1) 

RETURN (6-k_ + l)-(12-k_ + l)-(18-k_ + 1) 
k_ :+ 1 

f 10 10 

#28: Lp := fi nd_CN(RANDOM(10 )), q := fi nd_CN(RANDOM(10 )), GCD(p, q) , n := p-q, 

carmi chael_X(n) J 

#29: [p := 1299963601, q := 9624742921, 1, n := 12511815466282418521, 140400] 

T 16 

#30: Le := 2 + 1, d := INVERSE_MOD(e , 140400), dp := M0D(d , p - 1) , 

dq := M0D(d , q - 1) ] 

#31: [e := 65537, d := 68273, dp := 68273, dq := 68273] 

#32: [DIM(encrypt(t)) , timing(RANDOM(0) , encrypt(t))] = [190, 0.08s] 

#33: SOLVE(decrypt(encrypt(t)) = t) = true 

#34: timing(RANDOM(0) , decrypt(encrypt(t))) = 0.44s 

#35: FACTOR(n) = 601-1171-1201-1801.2341-3511 

We have already mentioned several times that anyone who can factor n, can also com- 
pute the secrret exponent d and hence decrypt any messages that use this RSA- 
scheme. Interestingly enough, the converse is also true: Anyone who knows d, can also 
factor n. In particular, two users should never share the same modulus n. 

The idea of the so-called common modulus attack is very simple. Take any x with 0<x<n 
and check if gcd(x,n)=l. If the unlikely case that x and n aren't coprime, this gcd(x,n) is 
one of the prime divisors of n and we are finished. Otherwise we have x A (ed) =x mod n 
and since x is invertible also x A (ed-l)=l mod n. Now looking at the congruence y A 2=l mod 
n, on the one hand it has the obvious solutions y= ±1 as well as two other solutions ±u, 
where u is a solution of the system of congruences x=l mod p, y=-l mod q using the CRT. 
On the other hand, we know that y=x A ((ed-l)/2) mod is also a solution. If y=±u mod n, 
this obviously implies that y±l is either divisible by p or q, but not both, and hence n can 
be factored into n=gcd(y+l,n)gcd(y-l,n). If y=-l mod n, we were unlucky and must choose 
another x. The same goes for the case, where y=l and (ed-l)/2 is odd. However, if y=l 
mod and (ed-l)/2 is still even, on can replace y by y=x A ((ed-l)/4) mod n and repeat as 
above. As you can see the chances of a success for each try are 50% and the number of 
tries is k, if 2 A k is the biggest power of 2 that divides ed-1. 
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Okay, let's try this out for a new standard setup of RSA. (The number in the textbox is 
the number of random values of x that were needed for a success. It is often 1 and 
rarely bigger than 2.) 

#36: setupO = true 

fi ndf actor (n , e, d, c_ := 0, k_, x_, y_) := 

Loop 

x_ RANDOM (n) 

y_ := GCD(x_, n) 

If y_ > 1 

[DISPLAY(c_) , RETURN y_] 
c_ :+ 1 

#37: k_ := e-d - 1 

Loop 
k_ :/ 2 

y_ := M0DS(x_ A k_, n) 

If y_ t ±1 

[DISPLAY(c_) , RETURN GCD(y_ - 1, n)] 

If y_ = -1 v (y_ = 1 a 0DD?(k_)) exit 


#38: fi ndfactor(n , e, d) = 

1813178269192617663103865449435130243024150198769632074073503976900611164996- 

3854618333374751591222286492608664190901804009784939978589851055114333847621- 

389 


#39: fi ndfactor(n , e, d) = 

6556616639787462592763064194321098871915177239834438525857041211204202630274- 

2297981875002873189458022619691366502304442132361484618947489970815071021403- 

29 

Okay, as we saw factoring n and finding the secret exponent d are computationally 
equivalent. On the other hand, we don't know whether the problem of RSA-decryption 
is as hard as the problem of factoring n (or finding d for that matter). There is a vari- 
ant of RSA, the so-called Rabin variant, where we can prove this though. For this vari- 
ant we use the exponent e=2, which is actually forbidden in RSA since the mapping x- 
>x~2 mod n is no longer a permutation of Zn. As a rule, if y is a square mod n, the equa- 
tion y-x^2 mod n has got 4 solutions. They can be found by finding first the solutions ±u 
mod p and ±v mod q, respectively, and then solving x=±u mod p, x = ±v mod q for all 4 
possible choices of signs, again by using the CRT. Usually only one of the 4 solutions 
makes sense when applying the source decoding, which enforces its uniqueness. 


Let's conclude with a demonstration of the Rabin variant using again an RSA-setup, set- 
ting e:=2 thereafter and simply ignoring the values referring to d. 


#40: [setupO, e := 2] = [true, e := 2] 


#41: c := encrypt(Let's call it a day!) 








